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SOLID  AEROSOL  PROGRAM 


1. 

INTRODUCTION 

1.1 

Future  NBC  (nuclear,  biological  and  chemical)  threats 
perceived  to  endanger  armed  forces  on  the  battlefield  may  Include  solid 
aerosols  as  well  as  the  more  familiar  vaporous  and  liquid  compounds. 

Protection  against  all  such  threats  Is  provided  most  effectively  by  collective 
protective  enclosures  supplied  with  filtered  air  at  a  suitable  overpressure, 
thus  preventing  contaminants  from  penetrating  against  the  pressure  gradient. 
The  effectiveness  of  enclosures,  whether  overpressure-protected  or  not.  In 
preventing  an  unsafe  level  of  exposure  is  clearly  dependent  on  controlling 
leakage  areas  through  which  contaminants  may  enter. 

Exposure  Is  the  integral  of  the  threat  concentration  over 
time.  The  measure  of  protection  provided  Is  defined  by  the  protection  factor, 
being  the  ratio  of  the  exposure  one  would  sustain  outside  the  enclosure  to 
that  sustained  Inside  the  enclosure  for  a  given  time  duration.  By  definition, 
the  protection  factor  le  greater  than  or  equal  to  one. 

The  determination  of  protection  factors  for  enclosures 
such  as  vehicles  and  shelters  is  based  to  date  on  modell.ing  the  behavior  of 
chemical  agent  simulant  vapors.  The  question  asked  here  Is  what  changes  may 
be  needed  to  protection  factor  models  based  on  vapors  to  account  for  size- 
dependent  mechanisms  affecting  the  defense  against  very  small  solid  particles? 
Further,  If  vehicles  and  shelters  need  to  be  tested  specifically  for  solid 
aerosols,  what  might  be  appropriate  integrity-testing  techniques? 

The  solid  aerosols  could  conceivably  take  any  of  the  following  forms j 

-  all  possible  NBC  agent  particles  (the  most  threatening  in  terms 
of  respiratory  penetration  being  in  the  range  0.5  -  5.0  microns) 
smoke,  dust,  obscurants,  and  resuspended  particulates 

The  threatened  objects  as  defined  for  the  research  program  could  involve i 

combat  and  armored  vehicles,  static  or  moving,  with  or  without 
cannon  fire 

-  fixed  or  mobile  shelter  structures 

-  either  leaky  or  over-pressure  protected 
suspect  filter  and  housing  Integrity 

-  various  leak  characteristics  (  smooth,  rough,  sticky,  greasy,  etc.) 

The  goals  of  this  research  program  are  to  develop  the 
methodology  needed  to  determine  solid  aerosol  protection  factors  for 
enclosures  in  the  particle  size  ranges  of  interest  and  to  express  the 
controlling  mechanisms  in  model  calculations.  At  a  later  stage  of  the 
program,  these  would  be  used  in  developing  integrity-testing  techniques. 


2 


1.2  Program  Outline 

1.2.1  Objectives 

The  research  reported  in  this  document  concerns  an  exploratory  program  executed  in 
the  framework  of  the  main  author's  one-year  term  as  an  Exchange  Scientist  with  the  Ventilation 
Kinetics  group  at  ERDEC  (Edgewood  Research ,  Development  and  Engineering  Center).  Given 
the  time  limitation,  on  the  one  hand,  and  the  generality  of  the  program  goals,  on  the  other,  there 
was  a  need  to  define  objectives  likely  to  be  achieved  under  the  constraints  of  time  and  knowledge 
gaps,  The  objectives  which  formed  the  basis  for  the  one-year  experimental  program  are  to: 

-  identify,  acquire  and  test  instruments  for  disseminating  and  sampling  candidate  solid 
aerosols 

-  generate  a  theoretical  understanding  of  physical  phenomena  to  enable  analysis  of 
experiments 

-  interact  with  other  individuals  and  agencies  to  learn  the  state  of  the  art  and  to 
prevent  duplication  of  efforts 

-  construct  an  experimental  facility  and  carry  out  trials  on  enclosures  to  determine 
vapor  and  solid  aerosol  protection  factors 

-  verify  the  results  in  light  of  current  theoretical  understanding  through  model 
calculations 


1.2.2  Experimental  and  Computational  Programs 

The  experimental  program  was  designed  to  exploit  as  much  of  the  existing 
capabilities  of  the  Ventilation  Kinetics  group  as  possible  while  keeping  in  mind  time  and  logistic 
constraints.  Furthermore,  after  clarifying  the  status  in  the  current  literature  on  enclosure 
protection  from  solid  aerosols,  the  decision  was  taken  to  concentrate  on  developing  the  ability  to 
measure  the  relevant  parameters  individually.  This  would  lead  to  a  better  fundamental 
understanding  of  the  controlling  phenomena  than  would  a  program  of,  say,  testing  integral 
exposures  of  actual  vehicles  or  shelters.  This  would  provide  the  foundation  for  future  stages  of 
the  program  once  areas  both  of  uncertainty  and  of  special  concern  were  identified  and  time  and 
logistic  constraints  could  be  lifted.  This  being  the  situation,  an  experimental  program  was 
constructed  with  the  following  guidelines; 

-  the  object  of  exposure  experiments  will  be  a  'dummy'  enclosure  with  a  volume  on 
the  order  of  a  typical  vehicle  cabin;  maximum  control  over  residual  leakage  areas 
should  be  ensured 

-  a  dedicated  test  facility  should  be  located  close  to  the  laboratories  of  Bldg.  3161 

-  for  greatest  flexibility  the  powder  dissemination  system  should  be  designed  to 
produce  a  cloud  of  infinitely-variable  duration  in  a  confined  but  open-ended  volume 

-  instrumentation  must  provide  real-time  measurement  and  computerized  data 
acquisition  for  all  relevant  physical  parameters 

-  candidate  aerosols  are  to  be  environmentally  acceptable 


3 


The  program  makes  use  of  accumulated  experience  in  the  application  of  gaseous 
tracers,  such  as  SF6,  in  the  determination  of  enclosure  air  exchange  late  and  as  a  vapor  challenge 
simulant.  SFg  is  an  internationally  recognized  tracer  for  air  movement  and  an  ideal  tool  for 
measuring  the  protection  factors  for  challenge  vapors.  This  technology  provides  the  basis  for 
incorporation  of  the  solid  aerosol  mechanisms,  since  air  movement  is  the  carrier  of  the  particulate 
threat. 


The  experimental  program  comprised  measurements  of  external  and  '  mal  size- 
dependent  solid  aerosol  count  rates  and  mass  concentrations  and  simulant  vapor  con;  .ntrations  as 
well  as  the  physical  parameters  needed  to  characterize  the  ventilation  kinetics  of  the  enclosure. 
Penetration  into  the  completely  sealed  enclosure  through  a  well-defined  aperture  was  driven  by 
pressure  differential.  Measurements  to  date  sufficed  to  permit  calibration  and  verification  of  a 
theoretical  model  of  vapor  and  solid  aerosol  protection  factors  as  envisioned  by  the  goals. 

Computations  and  models  developed  in  the  research  program  are  performed  in 
several  scientific  languages  for  IBM-compatible  personal  computers  acquired  especially  for  this 
purpose: 


- Mathcad  3.1  (Mathsoft,  Inc.) mathematical  modelling 

-  ASYSTANT  (Asyst  Software  Technologies,  Inc.)  numerical  calculations  and 
model  equation  curve  fitting 

-  TableCurve  3.10  (Jandel  Scientific)  empirical  data  curve  fitting 

-  Easyest  LX  (Asyst  Software  Technologies,  Inc.)  data  acquisition 

-  Graftool  (3-D  Visions  Corp.)  graphical  analysis 

Several  programs  were  written,  primarily  in  Mathcad  to  serve  the  research  program 
as  utility  codes  for  repetitive  operations.  These  include: 

-  background-corrected  particle  count  spectra 

-  concentration  decay  rates  (for  particles  and  vapors) 

-  theoretical  size-dependent  particulate  settling  rates  (calculation  and  optimization) 

-  optical  sensor  size-dependent  response  calibration 

-  size-dependent  particle  penetration  ratios 

The  final  product  of  the  computational  research  program  is  an  interactive  Mathcad 
model  for  the  calculation  of  the  vapor  and  size-dependent  solid  aerosol  protection  factors  for  a 
given  enclosure  and  aperture.  Penetration  dynamics  in  the  model  arc  driven  by  pressure 
differential  across  the  well-defined  aperture.  Measurements  of  the  ventilation  kinetic  parameters 
in  the  experimental  system  are  used  to  verify  the  air  exchange  rate  calculated  by  the  model. 
Theoretical  particle  ftltration  and  settling  rates  are  calculated  based  on  models  published  in  the 
literature.  Empirical  settling  rates  obtained  from  the  experimental  setup  can  also  be  used  to  verify 
measured  protection  factors.  The  model  is  a  flexible  tool  for  the  analysis  of  empirical  results, 
parameter  sensitivity  studies,  testing  procedure  design,  and  protection  estimation.  Results  of  the 
model  to  date  correlate  reasonably  well  with  measured  data  when  empirical  settling  rates  are  used. 
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1.2.3  Experimental  System  Requirements 

The  solid  aerosol  dissemination  system  is  a  critical  component  of  the  experimental 
research  system.  It  is  required  to  continuously  disseminate  a  variety  of  dry  powders  in  a  size 
range  on  the  order  of  0. 1  -  10  microns.  Energetic  de-agglomeration  is  necessary  to  ensure 
reproducible  size  spectra  for  the  smallest  of  sizes  with  as  little  pre-dissemination  treatment  as 
possible.  Mass  concentrations  in  the  range  1-100  mg/m^  are  required  to  be  maintained  outside 
the  test  enclosure. 

The  instrumental  setup  is  required  to  provide  data  on  the  following  parameters; 

-  air  exchange  rate  (  SFg  tracer  decay  method  with  one  minute  resolution) 

-  pressure  differential  (to  1  mm  Water  Gauge) 

-  air  flow  rate  (10  -  lOOO  1pm) 

-  temperature  differentials  (0.5  degree  C  resolution) 

-  external  and  internal  particle  mass  concentrations  (at  least  10^  inside/outside  ratio 
and  better  than  2  second  response  time) 

-  external  and  internal  particle  size  counts  (less  than  five  minutes  response  time  for 
one  spectrum) 

-  external  and  internal  SFg  tracer  vapor  challenge  concentrations  (at  least  10^ 
inside/outside  ratio) 

-  environmental  conditions  (wind  speed,  humidity,  etc.) 

1.3  Technical  Report  Objectives 

This  report  is  a  compilation  of  progress  achieved  during  the  research  program 
carried  out  from  March  1992  to  March  1993  at  ERDEC,  Starting  with  the  broadly  defined  goals 
as  stated  in  the  Introduction,  the  research  proceeded  systematically  in  an  effort  to  provide  a  sound 
technical  base  with  which  to  address  specific  questions  in  the  future.  It  has  been  a  project  of 
exploration  and  integration,  the  open  literature  concerning  enclosure  protection  from  solid 
aerosols  having  been  found  to  leave  many  unanswered  questions. 

Since  only  limited  practical  feedback  from  potential  users  of  the  project  at  ERDEC 
has  been  forthcoming,  the  results  have  remained  quite  general.  However,  to  the  extent  that  the 
experimental  facility  and  the  computational  tools  which  are  the  final  product  of  the  project  are  to 
be  exploited,  the  Technical  Report  will  aid  the  next  generation  of  researchers  in  absorbing  and 
employing  them.  The  scientific  background,  gleaned  from  the  literature,  is  included.  Each 
segment  of  the  experimental  setup  is  described  and  the  nuances  of  its  application,  as  experienced 
by  the  main  author,  is  explained.  Computer  programs  of  general  value  are  documented  for  future 
applications.  Their  use  requires  familiarity  primarily  with  the  Mathcad  3.1  scientific  programming 
language,  whose  use  is  highly  recommended. 

It  is  expected  that  the  fruits  of  this  program  will  aid  threat  analysts,  operational 
scientists,  simulation  and  modelling  scientists  in  addressing  their  questions  concerning  the  solid 
aerosol  challenge  with  greater  focus  than  previously. 


\\ 
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2.  LITERATURE  SURVEY 

2.1  \:aDQr  Challenge  Enclosure  Protection 

The  contamination  of  a  leaky'  enclosure  by  an  external  toxic  vapor  cloud  has  been 
treated  extensively  in  the  literature  (for  example;  Birenzvige,  1983a,b;  Wadden  and  Sheff,  1983; 
Chester,  1988;  Lewis,  1991).  Typical  results  for  some  pollutants  in  indoor  air  quality  studies  have 
been  summarized  by  Engelmann  (1990).  In  general,  the  enclosure  is  treated  as  a  well-mixed  single 
compartment  system  exposed  to  a  constant,  spatially  uniform  threat  concentration.  Analytical 
solutions  of  the  contaminant  balance  for  the  simplest  cases  provide  usefiil  equations  for  estimating 
the  efficiency  of  protective  measures.  These  solutions  are  often  given  in  terms  of  an  inside  to 
outside  ratio  (a  dose  reduction  factor)  or,  interchangeably,  an  outside  to  inside  ratio  (a  protection 
factor). 


The  vapor  protection  factor,  for  an  enclosure  is  defined  as  the  ratio  of  the  dose 
(vapor  concentration-time  integral)  whJch  would  <esult  from  exposure  to  an  outside  concentration 
during  time  T  to  the  dose  accumulated  inside  up  to  some  time  t  >  T.  The  governing  parameter  is 
volumetric  air  exchange  or  infiltration  rate,  R.  The  air  infiltration  rate  of  a  leaky  enclosure  is 
determined  by  a  complex  combination  of  factors  involving  the  nature  of  the  enclosure  envelope 
and  its  interactions  with  the  contaminated  environment.  It  represents  the  rates  both  of  penetration 
of  vapor  and  its  flushing  out  after  the  cloud  has  passed.  The  estimation  of  air  infiltration  into 
enclosures  in  the  short  and  long  term  is  dealt  with  extensively  in  the  literature  of  energy 
conservation,  indoor  air  quality,  and  environmental  health  and  is  well  beyond  the  scope  of  this 
discussion. 


The  degree  of  mixing  within  an  enclosure  is  also  a  complex  function  of  what  in 
many  cases  may  be  imponderables  unless  well-focused  measurements  support  its  estimation 
(Shair,  1974).  The  existence  of  dead  volumes,  channeling,  and  short  circuits  may  well  negate  the 
assumption  of  ideal  mixing  Suffice  it  to  say  that  estimates  of  the  air  exchange  rate  and  degree  of 
mixing  of  the  enclosure  in  question  for  the  time  period  involved  in  the  analysis  of  vapor  protection 
are  critical.  Estimates  may  possibly  derive  from  tracer  gas  decay  measurement  or  direct  air 
flowrate  when  feasible. 

Our  concern  here  includes  scenarios  in  which  the  exposure  duration,  T,  may  be 
short  compared  to  the  time  during  which  occupants  are  exposed  to  the  threat  inside  the  enclosure. 
Hence  the  protection  factor  must  account  for  both  time  periods.  Most  analyses  in  the  literature 
have  been  concerned  with  long-term  exposures  of  poorly  protected  enclosures  and  ignore  the 
exposure  subsequent  to  passage  of  the  cloud.  This  scenario  toils  under  the  assumption  that 
occupants  are  instantly  aware  of  the  passage  of  the  threat  and  immediately  exit  into  fresh  air,  an 
unrealistic  supposition  at  best.  A  full  derivation  of  the  equation  for  the  vapor  protection  factor 
providing  the  required  exposure  scenario  flexibility  is  given  in  Appendix  A.  For  a  uniform  non- 
cendensing  vapor  and  ideal  mixing  in  a  single  compartment,  the  balance  equation  leads  to; 
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i?- r+ exp(-i?  •/)•(! -exp(/2‘ 7^) 


Crack  filtering  and  interior  surface  deposition  of  particles  involve  three  removal 
mechanisms:  diffusion,  gravitational  settling  and  turbulent  impaction,  If  a  size-dependent  filtering 
factor,  /p,  and  settling  rate,  Pp,  are  included  with  the  air  exchange  rate  as  additional  removal 

terms  in  the  single  compartment  balance  model,  a  complete  expression  for  the  particulate 
protection  fkctor,  PFp,  is  obtained  (see  Appendix  A  for  the  fUll  development); 


R  -T 


where: 


.  7+ exp(-/?^ ./).(!-  exp(i?^  •  T))) 

Rp  =  R  +  Pp 


A  search  of  the  open  literature  relevant  to  sheltering  from  particulate  matter 
revealed  few  comprehensive  treatments  of  the  problem.  The  work  performed  by  Engelmann 
(1992a,b)  concerning  sheltering  effectivenes.s  of  buildings  and  vehicles  defines  a  dose  reduction 
factor,  the  inverse  of  the  protection  factor,  in  terms  of  particulate  removal  phenomena.  Here 
equilibrium  dynamics  for  a  monodisperse  challenge  were  studied,  but  underlying  size-dependent 
mechanisms  were  not  investigated.  In  the  following  sections  published  works  are  discussed  which 
deal  with  these  mechanisms  separately  and  which  are  here  subsequently  integrated  into  a 
comprehensive  model  for  calculating  the  size-dependent  particulate  protection  factor  for  an 
enclosure  given  in  Equation  2. 


The  rate  of  air  exchange  in  a  single  well-mixed  compartment,  R  in  equations  (1)  and 
(2),  is  estimated  ideally  by  the  ratio  of  the  air  flowrate  through  the  enclosure  to  the  volume  of  the 
enclosure; 


Standard  methods  for  measuring  the  air  exchange  rate  when  the  infiltration  flow,  Q, 
is  not  measurable  are  given  in  ASHRAE,  (1981).  A  practical  method  is  by  tracer  gas  dilution, 
covered  in  ASTM  Standard  E  741-80.  Tracer  gas  concentration  response  analysis  is  also  useful 
for  detecting  the  numerous  aberrations  from  ideal  mixing  which  may  occur  (Nauman,  1981). 
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When  considering  a  well-defined  aperture  as  the  sole  source  of  infiltration  into  an 
enclosure,  the  crack  flow  equation  can  be  used  to  estimate  the  air  flow  (ASHRAE  1981).  With 
differential  air  pressure  across  the  aperture  as  the  driving  force  for  penetration,  the  air  flowrate  is: 

C  =  (5) 

where  dp  is  the  differential  pressure,  n  is  the  characteristic  flow  exponent  which  varies  between 
0.5  and  1 .0  depending  on  the  Reynold's  number  ( Kreith  and  Eisenstadt,  1957),  C  is  a  discharge 
coefficient  and  A^  is  the  'effective*  leakage  area.  It  should  be  kept  in  mind  that  the  effective 
leakage  area  is  defined  by  the  above  relationship  between  the  pressure  differential  and  the 
flowrate.  It  will  not  in  general  equal  the  geometric  area  of  the  aperture  opening,  especially  if  the 
aperture  is  thick,  narrow,  rough  or  any  combination  of  these  factors. 

2.4  Stirred  Settling  of  Particles  in  Enclosures 

Solid  aerosol  settling  in  enclosures  with  convective  stirring  has  been  researched 
widely,  with  several  efforts  at  verifying  model  equations  through  controlled  chamber  tests  seen  in 
the  open  literature.  Activity  in  this  area  has  appeared  in  many  fields:  nuclear  safety,  fine  powder 
processing,  aerosol  sampling,  indoor  air  quality,  and  tropospheric  photo  chemistry. 

Earlier  works  strove  to  give  theoretical  meaning  to  the  observed  exponential  decay  in  the  particle 
count  concentration  in  small  chambers; 


dC(t)/cll=^-PpC(t) 


(6) 


where  the  decay  rate  comprises  uncoupled  stirred  settling  (j3^)  and  difflisive  deposition  ( /3^) 
components  (Van  de  Vate,  1972;  Harrison,  1979).  Their  analysts  led  to  the  expression: 


V  SD 
h  V'O 


(7) 


where  v  is  the  terminal  particle  settling  velocity,  h  is  the  chamber  height,  S  the  surface  area,  V  the 
volume,  D  the  particle  diflusivity,  and  a  the  difilision  boundary  layer  thickness  at  the  wall. 

Crump  and  Seinfeld  (1981)  later  derived  a  more  general  formula  for  the  turbulent 
deposition  and  gravitational  settling  of  particulates  in  arbitrarily  shaped  chambers  which  couples 
the  two  components  of  the  settling  coefficient  to  the  shape  of  the  chamber  and  gives  explicit 
expression  to  turbulent  eddy  diflusivity  for  a  stirred  chamber,  similar  to  the  results  of  Comer  and 
Pendlebury  (1951). 

Crump  and  Seinfeld's  equation  for  a  rectangular  chamber,  of  particular  interest  for 
our  project,  is,  as  expressed  by  Chen  et  al.  (1992a)  : 
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where  1,  w  and  h  are  the  length,  width  and  height,  respectively,  of  the  enclosure;  is  the  particle 

settling  velocity,  and; 

n-U, 

AT  - - 1 - —  (9) 

where  Dg  is  the  Brownian  difRjsion  coefficient  and  and  n  are  related  to  the  turbulent  energy 
dissipation  rate  (and,  hence,  to  the  rate  of  power  consumption/unit  mass  of  air)  due  to  all  sources. 
The  turbulent  eddy  diffiision  c'^efficient,  which  does  not  appear  explicitly  in  the  model,  has  been 
defined  as: 


De  =  ke>Y»  (10) 

where  Y  is  the  distance  fhom  the  wall.  The  model  of  Crump  and  Seinfeld  has  been  adopted  here 
for  use  in  the  model  calculations  of  settling  in  the  stirred  enclosure. 

It  v^ll  be  seen  later  that  the  two  parameters,  kg  and  n,  are  critical  to  the 
determination  of  the  particulate  settling  rates.  Several  attempts  to  quantify  and  correlate  them 
with  empirical  data  were  seen  in  the  literature.  Okuyama  el  al.  (1986)  reported  experimental 
observations  of  particle  settling  in  stirred  tanks  under  various  conditions  of  turbulent  energy  input. 
Results  were  compared  to  the  theoretical  results  of  Crump  and  Seinfeld  in  terms  of  particle  size 
and  turbulent  intensity.  The  exponent  n  has  a  theoretical  value  of  2  according  to  Crump  and 
Seinfeld,  but  Okuyama's  data  indicated  a  value  of  2.7. 

Chen,  et  al  (1992a)  used  solid  particle  decay  data  in  an  environmental  reaction 
chamber  to  determine  both  k^  and  n  as  a  function  of  vertical  temperature  gradients.  Although  the 
data  of  other  researchers,  which  Chen,  et  al  (1992a)  analyzed,  yielded  best  fits  for  n  in  the  range 
1.8S  -  3.S,  his  results  showed  n’»2  to  be  the  best  value,  as  Crump  and  Seinfeld  proposed. 
However,  the  wide  discrepancies  in  the  published  experimental  values  of  the  k^-n  pair  of 
parameters  sheds  considerable  doubt  on  the  generality  of  their  conclusion. 

2.5  Crack  Filtration  of  Solid  Aerosols 

The  literature  provides  little  comprehensive  data  on  the  size-dependent  infiltration  of 
solid  aerosols  into  structures  although  several  efforts  have  been  made  in  the  past  to  quantify  this 
factor  (Megaw,  1962;  Cristy  and  Chester,  1981;  Roed  and  Cannell,  1987).  For  example,  based 
on  early  results  Engelmann  (1990),  when  calculating  dose  reduction  factors,  assumed  that 
particles  below  2  microns  are  not  filtered  during  penetration  into  buildings.  Much  work  is  being 
in  recent  years  on  penetration  of  fine  particles  through  orifices  (Sutter,  et  al,  1980;  Mitchell,  et  al, 
1990;  Fletcher  and  Verkouteren,  1992),  although  under  larger  pressure  differences  than  interest 
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us  here,  and  through  respirator  facepieces  (Hinds  and  Kraske,  1987;  Hinds  and  Beilin,  1987; 
Holton  et  al,  1987;  Chen,  et  al,  1992b). 

A  review  report  by  Schwendiman  and  Sutter  (1977),  concerning  the  transport  of 
particles  through  gas  leaks,  summarized  equations  for  the  three  controlling  mechanisms: 
deposition  by  molecular  diffusion,  impaction  by  turbulent  eddies  and  gravitational  settling. 
Deposition  by  molecular  difhision  is  controlled  by  the  diffusion  coefficient,  aperture  length  and 
the  air  flow  rate.  Empirical  equations,  verified  for  diffusion  batteries  by  Sinclair,  et  al  (1976), 
were  presented.  Eddy  transport  of  particles  to  the  aperture  walls  is  determined  by  the  inerti^ 
properties  of  the  particles  in  the  turbulent  field.  The  fraction  lost  to  turbulent  impaction,  Fj,  is 
defined  by; 


(11) 

where  P  is  the  aperture  perimeter,  A  is  the  cross  sectional  area,  1^  is  the  length  of  the  aperture 
path,  V  is  the  average  velocity  and  K  is  the  true  deposition  velocity.  Lacking  sufficient 
experimental  data  for  the  very  narrow  apertures  of  interest  to  Schwendiman  and  Sutter  (and  to 
us)  required  to  correlate  K  with  particle  size,  Reynold's  number  and  surface  roughness,  the 
terminal  settling  velocity  is  used,  resulting  in  a  minimum  wall  loss. 

An  approximation  for  the  fraction  lost  by  gravitational  settling  in  a  rectangular 
aperture,  stipulating  a  uniform  velocity  profile,  is  given  by; 

Fo  =  ^  (IJ) 


2.6  Survey  Conclusions 

The  nature  of  the  research  program  was  determined  after  reviewing  the  'state  of  the 
art*  regarding  protection  of  enclosures  from  solid  aerosol  threats.  Given  the  absence  of  any 
comprehensive  treatment  of  a  size-dependent  particulate  protection  factor,  the  decision  was  taken 
to  integrate  available  knowledge  of  the  relevant  physical  mechanisms  into  a  flexible  mathematical 
model  which  would  become  an  analytical  tool  to  serve  the  program  goals.  In  parallel  with  the 
modeling  effort,  the  experimental  system  would  be  designed  to  provide  real  data  on  all  the 
parameters  required  to  calculate  the  protection  factors  for  both  vapors  and  solid  aerosols.  Data 
on  actual  protection  factors  would  serve  to  verify  the  model  for  the  exposure  scenarios  of  greatest 
interest. 


This  approach  was  chosen  over  other  options,  such  as  direct  measurements  of 
integral  vehicle  exposures  or  equilibrium  level  exposure  of  enclosures  (see  Engelmann,  1992b), 
not  just  for  short  term  logistic  reasons,  but  primarily  to  provide  a  tool  which  would  advance  the 
fundamental  understanding  of  physical  mechanisms  which  combine  to  determine  the  protection 
factor  for  any  scenario.  It  was  felt  that  direct  application  of  published  parameter  values  without  a 
firm  scientific  basis  could  yield  'blind'  results  of  limited  validity  and,  hence,  worth. 
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The  main  equations  for  the  protection  factor  model  as  applied  to  vapor  and  particles, 
with  variable  threat  duration  and  exposure  time,  are  shown  in  Appendix  A  as  derived  by  the  main 
author,  They  may  be  fed  either  theoretically  derived  parameter  values  or,  alternatively,  any 
combination  of  empirically  obtained  values.  This  includes  as  well  either  calculating  the  air 
exchange  rate  fhom  the  pressure-driven  flow  through  the  aperture  or  measuring  the  air  exchange 
rate  by  tracer  gas  decay. 

The  model  equation  of  Cttimp  and  Seinfeld  for  particulate  decay  in  a  rectangular 
enclosure  as  given  by  Chen  et  al.  (1992a)  has  been  adopted  in  the  overall  model  program,  The 
model  has  both  ke  and  n  as  independent  parameters.  In  view  of  the  wide  discrepancies  seen  in  the 
literature  regarding  these  parameters,  it  was  clear  that  the  investment  of  considerable  effort  in 
obtaining  specific  experimental  values  of  the  decay  rates  is  certainly  justified. 

Schwendiman  and  Sutter's  (1977)  report  served  as  the  basis  for  a  sub-model  to 
calculate  the  size-dependent  total  transport  fraction,  fp  in  Equation  (2),  of  particles  through  a 
well-defined  aperture.  Using  this  approach  the  theoretical  results  could  be  compared  to 
measurements  on  a  real  aperture  and  provide  viable  verification  of  the  theory  before  applying  It  to 
broader  questions.  Again,  this  is  as  opposed  to  the  option  of  attempting  to  study  penetration  on  a 
complex  system  such  as  a  real  vehicle.  Although  the  equations  are  relatively  simple,  ignoring  such 
possible  complications  as  electrostatic  charges,  aperture  clogging,  surfkce  roughness,  and 
tortuosity,  they  do  appear  robust  enough  to  provide  first-order  approximations. 


3,  DISSEMINATION  FACILITY 

3.1  Cloche  Construction 

Since  the  original  goal  of  the  program  is  to  investigate  solid  aerosol  penetration  into 
vehicles,  some  kind  of  field  testing  facility  was  envisioned.  Since  the  research  would  be 
exploratory  in  nature  and  difficult  to  execute  on  a  scheduled  basis,  the  facility  would  preferably  be 
independent  and  dedicated  to  the  program.  Logistically  speaking,  its  proximity  to  Bldg.  3161 
would  be  preferable,  allowing  use  of  its  laboratories.  An  enclosure  of  sufficient  dimensions  to 
house  vehicles  would  aid  in  confining  a  particulate  cloud  around  them  and  ensuring  uniformity 
and  relative  independence  from  the  vagaries  of  wind  conditions.  One  possibiiity  that  was  explored 
initially  is  construction  of  a  maintenance  tent  such  as  those  used  by  the  Vehicle  Recovery  ^^nit  on 
post.  Such  a  tent  has  been  ordered  and  eventually  could  be  raised  at  the  grenade  range  or  even  at 
M-Field  at  a  later  date 

Since  this  particular  solution  would  not  come  to  be  in  the  short  time  required  by  the 
program,  an  open-ended  agricultural  cloche  structure  identical  to  one  existing  behind  Bldg.  3161 
has  been  erected,  which  permits  controlled  dissemination  of  solid  aerosols  and  simulant  gases  and 
exposure  of  small  enclosures.  The  cloche  is  made  of  semi-circular  metal  arches  2  mt^tprs  high, 

4,3  meters  wide,  covered  with  polyethylene  sheeting  to  a  length  of  1 8  meters.  Plasticized  nylon 
tarps  provide  shade  over  the  work  areas  at  either  end.  A  black  sun-screen  netting  covers  the  rest 
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of  the  cloche  between  the  ends.  The  cloche  Is  anchored  on  the  concrete  pad 
area  adjacent  to  the  southern  end  of  bldg.  3161.  The  cloche  is  intended  for 
use  on  soil  surfaces  by  simply  driving  2“  metal  pipes  into  the  ground  to 
support  the  arches.  One  side  was  constructed  as  intended,  along  the  edge  of 
the  concrete  pad.  The  other  side  had  to  be  anchored  onto  the  concrete.  For 
this  purpose,  specially-made  metal  plates  with  sections  of  the  original  pipes 
welded  to  them  were  produced  and  anchored  to  the  pad  with  3/4"  bolts.  A 
photograph  of  the  completed  cloche  facility  is  provided  in  the  Appendix. 

The  dimensions  of  the  cloche  were  intended  to  permit 
dissemination  of  a  solid  aerosol  cloud  at  the  western  end  and  exposure  of  a 
chamber  or  small  vehicle  (a  HMMWV  vehicle  can  enter  the  cloche)  at  the  eastern 
end.  It  was  expected  that  the  18  meter  length  would  suffice  to  guarantee 
development  of  uniform  concentrations  in  the  vicinity  of  the  target  enclosure 
with  a  properly  designed  dissemination  system. 

The  facility  receives  electric  power  from  an  external 
receptacle  and  compressed  air  from  the  laboratory  compressor  through  a  1/2" 
plastic  pipe  leading  from  the  air  supply  receptacle  on  a  nearby  room.  Metal 
link  connectors  are  attached  to  the  plasticized  tarp  at  the  western  end 
permitting  an  additional  tarp  to  be  hung  there  as  a  rain  shield. 

The  facility  is  intended  to  provide  for  dissemination  of 
particulate  clouds  of  varying  duration,  from  short  pulses  of  about  one-half 
minute  to  extended  exposures  of  even  one  half  hour.  Hence  it  is  operated  with 
both  ends  open  so  that  the  dissemination  fan  system  also  flushes  the  facility 
after  the  passage  of  the  cloud.  If  so  desired  in  the  future,  the  ends  can  be 
closed  by  tarps  to  trap  the  cloud  and  produce  extended  stable  exposure  times. 

3.2  Solid  Aerosol  Dissemination  System 

The  dry  powder  dissemination  system  is  reauired  to  produce 
a  controlled,  continuous  stream  of  particles  which  have  been  de-agglomerated 
to  the  fullest  extent  possible.  It  would  potentially  be  used  for  primary 
particles  in  the  range  0.  02  -  30  microns,  hence  utmost  flexibility  in 
operating  conditions  is  necessary. 

Discussions  with  Mr.  R.  Doherty  of  the  Obscuration 
Sciences  Branch  led  to  the  suggestion  to  use  a  similar  component  to  that  which 
produces  particle  clouds  for  that  Branch's  aerosol  chamber  facility  -  the  SRI 
sonic  velocity  nozzle  (hereafter  referred  to  as  the  sonic  nozzle).  The 
component  is  simple  and  effective  as  a  de-agglomerating  apparatus.  Its 
application  is  well  documented  in  the  literature  (Withal  and  Gates,  1983).  A 
spare  unit,  whose  schematic  diagram  is  provided  in  Figure  1,  was  adopted  for 
use  in  the  dissemination  facility. 

The  sonic  nozzle  uses  air  flowing  at  the  speed  of  sound 
to  produce  a  region  of  high  shear  force  where  the  powder  enters  the  air 
stream.  The  flow  of  air  past  the  opening  of  the  powder  feed  tube  causes  a 
suction  which  draws  particles  into  that  region  where  the  shear  forces  and 
pressures  caused  by  fluid  drag  then  act  to  break  apart  agglomerates  comprised 
of  the  primary  particles.  Van  der  Waals  forces  are  the  primary  adhesive 
mechanism  forming  agglomerates  in  dry  powders.  The  energy  available  at  the 
point  of  pressure  drop  outside  the  nozzle  divides  between 
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suction  and  de-agglomeration  of  the  powder  drawn  into  the  region.  Operation  of  the  nozzle 
requires  a  compromise  between  the  two  mechanisms  to  suit  the  needs  of  the  application. 

At  this  point  it  remained  to  obtain  some  apparatus  for  controlled  continuous  feeding 
of  powder  to  the  sonic  nozzle  feed  tube.  The  solution  was  found  in  a  Metronics  Associates,  Inc. 
Fluorescent  Particle  Aerosol  Generator  (Model  10),  also  provided  by  the  Obscuration  Sciences 
Branch.  This  apparatus  is  intended  to  provide  a  continuous  release  of  dry  powder  particles  into 
the  atmospheric  environment  as  simulants  for  the  transport  and  dispersion  of  airborne  particulate 
contaminants.  Powder  is  fed  fbom  a  stirred  hopper  by  a  rotating  toothed  wheel  into  the  intake  of 
a  high  speed  centrifugal  blower.  The  feed  rate  is  controlled  by  the  rotational  speed  of  the  wheel. 
Two  alternative  wheels  with  different  tooth  spacings  are  available,  providing  two  broad  ranges  of 
feed  rates.  The  hopper  has  about  a  2  liter  capacity.  Since  the  toothed  wheel  feeds  powder  by 
discrete  increments,  the  mechanism  is  not  truly  continuous  but  pulsed.  This  was  not  expected  to 
be  of  concern  even  at  low  feed  rates  since  the  powder  stream  would  be  fUrther  dispersed  by  the 
turbulent  means. 

The  final  dissemination  system  is  then  comprised  of  the  aerosol  feeder  coupled  to 
the  pneumatic  nozzle,  the  blower  having  been  removed  (see  Figure  2  and  the  photograph  in 
Appendix).  The  nozzle's  powder  feed  tube  is  connected  to  a  smooth  tygon  tube  to  a  funnel 
attached  to  the  underside  of  the  powder  feeder.  The  nozzle  is  fed  by  the  laboratory  air 
compressor  operating  at  70  psi  with  a  3  psi  differential  range.  About  20  inWG  suction  develops 
at  the  powder  input  so  that  as  each  parcel  of  powder  drops  from  a  gear  tooth  into  the  funnel  it  is 
drawn  rapidly  into  the  nozzle  and  disseminated.  The  air  flow  ratio  between  nozzle  output  and 
powder  input  is  about  50:1,  so  the  shear  force  de-agglomerating  mechanism  is  considerably 
energetic. 


The  powder  feed  apparatus  is  located  on  a  table  at  the  head  of  the  cloche  together 
with  its  control  console.  The  nozzle  is  situated  on  a  small  table  in  front  of  the  feeder  at  1.25 
meters  above  the  ground.  Behind  it  are  three  24"  fans  in  a  row  facing  down  the  length  of  the 
cloche.  A  12"  fan  situated  in  front  of  the  nozzle  faces  into  the  stream  of  disseminating  particles, 
spreading  the  stream  across  the  plane  of  the  larger  fans  and  facilitating  its  uniform  crosswise 
mixing  (see  photographs  in  the  Appendix).  The  system  produces  a  1-1.5  m/s  breeze  at  a  distance 
of  1 S  meters  from  the  fans  near  the  far  end  of  the  cloche.  Particle  mass  concentration 
measurements  later  showed  that  deviations  of  only  10-20%  from  a  mean  value  across  the  area  of 
the  cloche  at  that  distance  may  be  achieved.  Cloud  concentrations  up  to  several  tens  of  mg/m^ 
can  be  produced,  depending  on  powder  composition,  feed  rate,  and  fan  power  level. 


Figure  2;  Dissemination  System  Schematic  Diagram 
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3.3  Tracer  Gas  Dissemination 

SFg  gas  is  used  universally  as  a  tracer  for  air  movement  and  the  transport  of  vapor 
and  gasses  in  general.  Its  use  as  a  means  for  measurement  of  air  exchange  rates  is  well 
documented  and  standardized  (ASHRAE,1981). 

The  gas  is  supplied  liquefied  under  a  vapor  pressure  of  about  20  atmospheres. 

Hence,  it  is  disseminated  under  its  own  pressure  through  a  regulator  and  flowmeter  of  some  sort. 
A  bottle  of  SFg  is  located  near  the  particle  dissemination  system  when  the  penetration  of  vapors  is 
to  be  investigated  with  the  SF5  simulant.  Dispensing  the  gas  through  a  rotameter  and  fine 
hypodermic  needle  as  a  controlling  orifice  has  proved  sufficient  for  our  purposes.  The  needle  is 
situated  near  the  outlet  of  the  sonic  nozzle  so  that  the  gas  experiences  the  same  amount  of 
turbulent  dispersion  as  do  the  particles. 


4.  TEST  CHAMBER 

4.1  Chambfir.Cflhslagtifln 

The  exposure  testing  enclosure  is  a  reconstructed  gas  mask  testing  chamber  made  of 
fiberglass  and  having  several  Plexiglas  windows  (see  photograph  in  Appendix).  The  airlock  wall 
was  left  out  to  make  a  volume  of  4.43  m^,  sufficient  to  repre.sent  a  small  crew  compartment.  The 
chamber  lies  on  its  side  on  a  small  wheeled  carriage  for  ease  in  movement.  All  joints  were  sealed 
with  1/4"  X  2"  rubber  stripping  to  ensure  air-tightness.  Provisions  were  made  for  several  signal 
cables  and  air  sampling  tubes. 

One  window,  facing  toward  the  dissemination  system,  has  been  replaced  by  two  10 
mm  thick  plates  of  perspex  held  to  the  rubber  seal  by  16  quick-release  clamps.  The  85  cm.  long 
plates  can  be  separated  using  spark  plug  feeler  gauges  to  provide  a  well-defined,  controlled  leak 
into  the  chamber.  Additional  perspex  strips  are  added  at  the  interface  to  provide  leak  thicknesses 
up  to  5  cm.  as  desired. 

4.2  Instrumentation  Description 

Figure  3  provides  a  schematic  outline  of  the  instrumentation  system  constructed  to 
collect  and  analyze  all  the  parameters  related  to  investigating  vapor  and  particulate  protection  of 
the  exposure  test  chamber.  The  chamber  is  instrumented  to  measure  inside  and  outside  particle 
mass  concentrations  and  size  count  and  vapor  concentrations  as  well  as  air  flow,  temperature  and 
pressure  differences,  and  relative  humidity. 

A  100  cfm  centrifugal  blower  at  the  far  end  of  the  cloche  is  connected  to  the 
chamber  through  3.5"  plastic  tubing  to  provide  over-pressure  or  under-pressure  as  needed.  The 
blower  output  is  connected  to  a  metal  iris  which  permits  continuous  adjustment  of  the  flowrate 
into  or  out  of  the  chamber.  A  Kurz  Model  SOS  mass  flowmeter  located  with  its  control  box  on 


Figure  3;  Solid  Aerosol  Instrumentation  System 
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the  flow  line  near  the  control  table  measures  air  flowrate.  The  pressure  differential  between  the 
chamber  and  the  environment  is  measured  by  an  Endevco  8510B  (1  psi  full  scale)  niezo-resistive 
transducer  located  with  its  amplifier  inside  the  chamber.  The  blower  produces  pressure 
differentials  of  several  mmWG  only.  If  higher  pressure  differentials  are  required,  a  400  cfrn 
blower  is  available,  but  structural  limitations  of  the  fiberglass  construction  and  the  rubber  strip 
sealing  of  the  windows  must  be  considered. 

Three  MIE,  Jnc.  RAS-2  Real-time  Aerosol  Sensors  situated  around  the  chamber 
obtain  cloud  mass  concentration  data  in  the  range  0.01  - 100  mg/m^  with  a  0.2  second  response 
time.  The  RAS-2  sensors  are  light-scattering  nephelometers  providing  analog  signals  linearly 
proportional  to  mass  concentration.  The  sensors  are  provided  with  both  open  (flee  convective) 
chambers  and  closed  (pumped)  chambers.  Experience  showed  that  the  closed  chamber  option 
maintained  the  sensor  windows  cleaner  than  in  the  open  chamber  option,  especially  when 
disseminating  talc  or  titanium  dioxide  powders  which  tend  to  stick  to  the  glass  windows  of  the 
sensors  in  humid  weather.  The  three  sensors'  chambers  are  connected  by  tygon  tubing  to  a 
manifold  and  a  10  1pm  diaphragm  pump  which  pulls  air  through  the  chambers.  Valves  on  each 
arm  of  the  manifold  allow  balancing  of  the  air  flow  from  each  of  the  sensors.  The  flow  is  adjusted 
to  about  2 1pm  tlu'ough  each  sample  chamber. 

Inside  the  enclosure  an  MIE,  Inc.  RAM-1  Real-time  Aerosol  Monitor  measures  the 
penetrating  particle  mass  concentration  in  the  range  0.001  -  200  mg/m^  with  response  time  down 
to  O.S  second.  This  monitor  uses  the  same  detector  system  as  the  RAM-2  sensor  but  provides 
digital  readout  of  the  concentration  as  well  as  analog  output,  internal  reference  calibration,  and 
continuous  flushing  of  optical  surfaces. 

A  Particle  Measuring  Systems  LAS-X  Laser  Aerosol  Spectrometer  inside  the 
enclosure  collects  particle  size  count  spectra.  The  LAS-X  provides  64  bins  in  four  ranges  in  the 
region  of  0.12  -  7.5  microns,  The  LAS-X  contains  a  precision  He-Ne  laser-illuminated  optical 
mirror  system  which  permits  single  particle  sampling  by  collecting  the  light  scattered  flom  each 
particle  which  is  then  sensed  by  a  solid  state  photodetector.  The  instrument  is  calibrated  by  the 
manufacturer  based  on  diameter  using  spherical  latex  particles.  The  data  output  is  sent  by 
standard  RS-232C  interface  signal  to  an  HP-8SA  portable  computer  for  processing  and  later 
analysis.  A  TSI  APS  3300B  particle  sizer  is  also  available  for  particle  sizes  up  to  30  microns  but 
has  not  yet  been  put  to  use  since  concern  to  date  has  been  with  the  0.1-5  micron  range. 

The  L.AS-X  spectrometer  input  is  connected  through  a  tygon  tube  to  the  outside  wall 
of  the  enclosure  when  the  external  cloud  spectrum  is  to  be  obtained.  In  this  configuration  the 
LAS-X  is  linked  to  a  TSI  Model  3302  Diluter  to  provide  lower  countrates  and  prevent  saturation 
in  the  lowest  channels.  The  maximum  countrate  in  channels  1-16  noted  by  the  manufacturer  is  7.0 
KHz/channel.  Since  the  sampling  flowrate  is  only  5-10  cc/sec,  use  of  the  TSI  3302  diluter 
necessitated  construction  of  a  special  flow-matching  arrangement  to  provide  isokinetic  sampling 
by  the  LAS-X.  This  apparatus  is  discussed  in  Appendix  E. 

The  data  provided  by  the  RS-232C  output  are  processed  in  the  HP-85A  computer  by 
a  dedicated  program  written  for  this  purpose  at  Los  Alamos  National  Laboratory  (Soderholm, 


18 


1984).  The  program  output  provided  on  the  HP-85  A's  printer  summarizes  information  relevant  to 
each  run.  These  include  a  plot  of  the  bin  count  data  for  each  of  the  four  ranges,  a  chart  of  the  run 
parameters,  and  a  listing  of  particle  counts  for  each  of  the  16  bins  in  all  the  four  size  ranges. 

Since  the  size  ranges  are  overlapped  and  not  of  equal  bin  width,  analysis  of  the  data  may  be 
confusing.  The  program  calculates  a  best-fit  log  normal  distribution  of  the  data  and  plots  it  with 
the  original  counts  on  the  output  graph.  This  calculation  is  useful  only  if  the  data  are  not  bimodal 
with  the  background  spectrum,  in  which  case  the  data  must  be  analyzed  separately  by  first 
subtracting  the  background  spectrum.  Manipulation  of  the  LAS-X  data  will  be  discussed  further 
in  Appendix  C.  1 .  A  sample  of  the  data  output  of  the  program  for  ambient  background  is  provided 
in  Figure  4.  The  data  plot  and  parameter  chart  are  shown. 

Several  options  exist  for  the  measurement  of  SF5  tracer  gas  concentrations  inside  and 
outside  the  enclosure.  A  MIRAN-1 A  Infrared  Gas  Analyzer  can  be  situated  inside  the  enclosure 
to  monitor  tracer  concentration  during  air  exchange  rate  measurements  of  the  sealed  system  or 
when  the  rate  is  relatively  slow.  Fast  response  measurements,  such  as  for  vapor  penetration 
testing  and  short  air  exchange  rates  must  be  done  using  an  electron  capture  gas  chromatograph 
and  automatic  sampling  system.  Such  a  system  has  been  acquired  from  HNU,  Inc.  but  has  yet  to 
be  adapted  to  the  overall  system.  Pending  its  inclusion  in  the  experimental  setup,  the  available  S- 
Cubed,  Inc.  (now  Lagus  Associates)  portable  SF5  Gas  Chromatographs  have  been  used.  The 
columns  have  been  changed  from  molecular  sieve  to  alumina,  permitting  more  rapid  analysis  of 
SFg.  The  SF5  peak  follows  the  oxygen  peak  and  clears  the  systems  in  less  than  one  minute.  For 
trial  runs  the  injection  cycle  has  been  hand  operated  from  the  instrument  panel.  A  10 1pm 
diaphragm  pump  at  the  control  table  samples  air  from  the  chamber  in  a  closed  cycle  while  the 
sampling  loop  of  the  gas  chromatograph  is  permitted  to  bleed  off  the  required  flow  during  the 
sample  injection  cycle.  The  vented  santple  air  is  returned  to  the  return  air  line  to  the  chamber, 
thus  preventing  even  the  slightest  effect  on  the  measured  air  exchange  rates. 
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DATE  VY/MM/DD  TIME  1€.351£  PSLF 
MAX  CNTS/SEi>  1£32  SEC»  8 
R>0-3  TRLIN  0  BINS  5.00  ML.-'SEi:: 
DILUTION  RATIO  «  1.000E+000 

TEMPERATURE  2.200E  +  001 

ATM  PRESSURE  <mm  H?>=  7  ,  i500E  +  002 
REL  HUMIDITY  <:i)=  5.000E-f001 

A  PARAMETER  =  0.000E+000 

B  parameter  =  0.000E+00O 

C  PARAMETER  =  0.000E+000 

HUM  CONC  <NUMB.-'cm3)=  2.510E  +  002 

GEOM  MEAN  DIAM  <:.um>=  1.746E-001 
GEOM  STANDARD  DEV  =  1 . 348E+000 
MASS  CONC  <m^/in3)=  1.879E-003 
UNDIL  MRSSCON<tn^- !n3)=  1.879E-003 
MAS  CHTSxSEC  =  I . 232E+003 

PEAK  DIAMETER  <vrO  =  1.230E-001 
FIT  HUM  CONC  <#.-'cm3)=  0.000E  +  000 
FIT  GEOM  MN  DIR  1  . 000E-008 

FIT  GEOM  STAND  DEU  =  1 , 000E+002 
FIT  UNDIL  M0  <m^/m3)=  0.000E+000 


Figure  4;  Sample  LAS-X  Data  Output  for  Ambient  Air 
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5 .  SOLID  AEROSOL  PROTECTION  FACTOR  MODEL 

The  protection  factor  model  program  is  written  in  the  Mathcad  3. 1  scientific 
programming  language.  This  is  an  interactive  tool  operated  in  the  Microsoft  Windows 
environment  which  permits  model  equations,  text,  and  graphs  to  be  written  in  the  same  program 
as  if  on  a  scratchpad.  As  described  earlier,  the  model  calculates  particle  size-dependent 
parameters  for  the  solid  aerosol  protection  factor  in  Equation  (2).  For  comparative  purposes,  it 
also  calculates  Equation  (1)  for  the  vapor  protection  factor.  The  program  is  also  provided 
experimental  data  for  particulate  settling  rates  and  calculates  a  semi-empirical  result  for 
comparison  to  the  fiilly  theoretical  results. 

After  presenting  the  size-dependent  protection  factors  in  a  two-dimensional  matrix 
for  several  pressure  differential  values,  the  results  can  be  weighed  either  by  a  given  experimental 
size  spectrum  or  by  a  log  normal  distribution  spectrum  defined  by  user-supplied  geometric  mean 
and  standard  deviation  parameters.  The  log  normal  weighted  results  can  be  displayed  with  both 
the  theoretical  settling  rates  and  the  experimental  settling  rates  alongside  the  vapor  protection 
factors. 


Since  the  program  links  text  and  model  equations  in  the  same  ftamework,  a  sample 
run  of  the  program  is  presented  in  Appendix  B  in  the  original  Mathcad  code  with  explanatory  text 
included. 


The  parameters  required  for  the  particulate  protection  factor  calculation  include 
enclosure  dimensions,  aperture  dimensions,  differential  pressure  (five  values  are  used 
simultaneously  to  provide  immediate  sensitivity  results),  turbulent  mixing  parameters  (ke  and  ne  in 
the  program),  challenge  and  exposure  durations  and  challenge  size  spectrum  within  a  0. 1-10 
micron  range.  The  parameters  used  for  the  sample  run  shown  in  the  Appendix  are  as  follows: 


-  enclosure  dimensions; 

-  aperture  dimensions: 

-  cloud  duration  (T): 

-  particle  density: 

-  log  normal  parameters 

-  flow  exponen  (n): 

-  turbulent  dissipation: 

-  pressure  differentials: 


1  =  213  cm.  h=117cm.  w=167cm. 

1  =  80  cm.  h=  0,015  cm.  w  =  4cm. 

5  minutes  exposure  time  (t):  30  minutes 

2,7g/cm3  (Talc) 

geometric  mean  1.2  Standard  deviation:  2.284 
0.5 

ke  =  1000  ne  =  2.038  (approx,  small  fan  operating) 
-1,  -1.5,  -2,  -2,5  and  -3  mmWG  (simulating  wind) 


These  parameters  were  chosen  to  demonstrate  a  situation  in  which  significant 
filtration  of  the  larger  particles  in  the  narrow  and  deep  aperture  should  be  assured.  The  low  air 
exchange  rates  expected  from  the  small  pressure  differentials  assure  sufficient  time  for 
pronounced  particulate  settling  to  occur. 


The  general  features  of  the  results  correlate  well  with  expected  behavior.  Particle 
filtration  in  the  aperture  (fig.  5)  increases  sharply  at  1  micron.  As  the  pressure  differential 
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increases  (higher  curves),  the  filtration  decreases  as  more  and  larger  particles  are  pulled  through 
the  aperture.  The  particle  settling  rates  (fig.  5)  show  a  minimum  at  the  junction  of  the  diffijsional 
and  gravitational  mechanisms  at  about  1.5  micron  in  this  case.  The  particulate  protection  factors 
given  in  Figure  7  along  with  the  vapor  protection  factors  (straight  lines)  show  pronounced  size* 
dependence  with  some  improvement  at  the  smallest  sizes  and  sharp  improvement  at  the  larger 
sizes.  This  behavior  is  a  combination  of  the  higher  decay  rates  at  the  smallest  sizes  and  high 
filtration  and  settling  rates  of  the  larger  sizes.  The  weighted  protection  factors  in  figure  6  for  the 
log  normally  distributed  spectrum  shown  in  Figure  8  are,  of  course,  correlated  with  the  vapor 
protection  factors  but  exhibit  higher  values  due  to  the  particulate  removal  mechanisms.  As  the 
differential  pressure  increases,  the  increased  air  exchange  rate  reduces  both  the  vapor  protection 
and  the  advantage  obtained  by  particle  filtration. 

The  same  functional  form  of  the  size*dependent  protection  factors  is  to  be  expected  in 
all  instances.  Their  values  can  only  be  larger  than  those  for  vapors  due  to  removal  mechanisms. 
The  extent  to  which  these  mechanisms  influence  the  protection  for  a  given  scenario  in  comparison 
to  the  vapor  protection  depends  on  the  solid  aerosol  size  spectrum.  The  example  shown  in 
Figures  9  relates  to  the  broadly  polydispersed  threat  shown  in  Figure  8.  Figures  10-15 
demonstrate  the  effect  on  the  weighted  protection  factors  of  hypothetical  relatively  monodisperse 
spectra.  The  spectrum  in  Figure  10  with  a  geometrical  mean  size  of  4  microns  yields  the 
weighted  factors  shown  in  Figure  1 1 .  Note  the  sharp  rise  compared  to  the  polydisperse  spectrum 
of  Figure  8.  The  influence  of  the  differential  pressure  is  also  accentuated  as  the  larger  particles 
are  more  strongly  filtered  at  the  lower  pressures.  The  1  micron  particles  of  Figure  12  track  the 
vapor  protection  factor  closely.  They  are  not  influenced  by  the  change  in  pressure  as  the  filtering 
factor  is  close  to  1.  Their  factors  are  the  closest  to  those  of  the  vapor  as  expected  by  their 
minimum  settling  rate.  The  smallest  particles,  shown  in  Figure  14,  show  the  effect  of  higher 
settling  rates  in  Figure  15. 

After  the  veracity  of  the  theoretical  model  is  shown  to  be  confirmed  through 
experimental  verification  in  the  coming  chapters,  further  sensitivity  analyses  will  demonstrate  the 
usefulness  of  the  model  program  as  a  tool  in  prediction  of  solid  aerosol  protection  and  in 
designing  of  integrity  tests. 
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Figure  5:  Total  Transport  Fraction 
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Figure  6;  Theoretical  Paiticle  Settling  Rates  [I /sec] 
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Figure  7:  Size-dependent  Particulate  and  Vapor  Protection  Factors 
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Figure  8:  Log  Normal  Challenge  Spectrum 


Figure  9;  Theoretical  Weighted  Protection 
Factors 
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Figure  10:  Log  Normal  Spectrum 
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Figure  12:  Log  Normal  Spectrum 
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Figure  14:  Log  Normal  Spectrum 
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Figure  11;  Theoretical  Weighted 
Protection  Factors 
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Figure  13:  Theoretical  Weighted 
Protection  Factors 
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Figure  15:  Theoretical  Weighted 
Protection  Factors 
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6.  EXJ'LtUMENTAL  MEASURE^ffii^^^S 

6.1  Einslaflicls  Spsclra 

Several  candidate  simulant  powders  were  examined  for  their  ease  in  dissemination 
and  usefulness  as  test  particles  within  the  0. 1-10  micron  range.  Pending  more  specific  criteria  in 
choosing  threat-specific  candidate  simulants,  only  powders  whose  use  in  the  open  air  is  acceptable 
by  environmental  safety  standards  were  considered.  Those  tested  include: 


-  Asbestos-^ee  Talc 

-  Aluminum  oxide 

-  Titanium  dioxides 

-  Hydrite  MP  clay 


-  Hydrophobic  Aerosils  (DeGussa  fUmed  silicas) 

-  Hydrophobic  Sipemats  (DeGussa  precipitated  silicas) 
-Coarse  Arizona  road  dust  (AC  dust) 

-Cabosils  (Cabot  fUmed  silicas) 


A  main  concern  is  the  efficiency  of  the  powder  feeder  and  sonic  nozzle  in  de- 
agglomerating  and  disseminating  bulk  candidate  powders.  Powders  have  been  tested  for  flow 
ability  in  the  feeder  system  and  their  disseminated  size  spectra.  As  discussed  above,  the  addition 
of  scraper  arms  on  the  stirring  rod  greatly  improved  the  efficiency  of  the  powder  feeder  apparatus 
except  for  the  most  agglomerated  powders  on  damp  days.  Example  size  spectra  for  bulk 
candidate  powders  following  dissemination  through  the  nozzle  and  transport  over  a  distance  of  1 5 
meters  are  shown  in  Figure  16.  These  spectra  were  obtained  fi'om  LAS-X  data  by  subtracting 
background  counts  in  a  dedicated  Mathcad  program  which  will  be  explained  fUrther  on  in  the 
Report.  It  was  hoped  that  a  wide  repertoire  of  solid  aerosol  spectra  could  be  constructed  with  as 
little  preliminary  powder  preparation  as  possible.  Both  narrowly-  and  broadly-dispersed  spectra 
should  be  available  in  a  wide  specific  mass  range. 


The  Aerosils,  Cabosils  and  Sipemats  have  primary  particle  sizes  in  the  7  -  20  nm 
range,  although  supplied  as  agglomerates  of  varying  sizes.  The  hydrophobic  Aerosils  are  basic 
fumed  silicas  which  have  been  treated  with  silane  to  produce  varying  levels  of  hydrophobicity, 
depending  on  the  primary  particle  size.  As  a  result  of  their  hydrophobicity,  they  are  very  readily 
disseminated  and  can  be  excellent  candidates  for  simulating  solid  aerosols  under  0.2  micron.  It 
should  be  kept  in  mind  that  the  lower  limit  of  detection  of  the  LAS-X  is  0. 12  micron,  so  that  it  is 
not  known  if  particles  are  produced  from  the  Aerosils  and  Cabosils. 


Certain  agglomerate  spectra  may  be  tailored  somewhat  by  changing  the  nozzle  energy 
or  by  bypassing  it  if  larger  sizes  are  needed.  In  that  case  the  original  blower  of  the  Metronics 
Aerosol  Generator  is  u  ^ed.  The  measure  of  de-agglomeration  achieved  in  the  nozzle  for  Aerosils 
and  Sipemats  is  shown  qualitatively  in  Figures  17  and  18.  In  Figure  17  the  Aerosil  812  spectra 
for  dissemination  through  the  original  blower  of  the  Metronics  system  is  shown  in  comparison  to 
that  obtained  fi'om  the  sonic  nozzle.  Essentially  all  the  agglomerates  above  0.2  micron  have  been 
separated,  by  the  nozzle.  In  Figure  18  Sipemat  D-1 1  spectra  are  shown  where,  in  addition,  the 
spectrum  obtained  by  simply  dispersing  the  bulk  powder  in  front  of  the  large  fans  is  also  shown. 
Agglomerates  up  to  4  microns  have  all  been  de-agglomerated  to  less  than  O.S  microns.  The 
Sipemat  D-1 1  is  supplied  as  a  nominal  4  micron  agglomerate. 
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Figure  16:  Fine  particle  spectra  from  sonic  nozzle  dissemination  (continued) 
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Figure  16:  Fine  particle  spectra  from  sonic  nozzle  dissemination  (continued) 
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Figure  16:  Fine  particle  spectra  from  sonic  nozzle  dissemination  (continued) 
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Figure  16;  Fine  particle  spectra  from  sonic  nozzle  dissemination  (continued) 
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Figare  16;  Fine  particle  spectra  from  sonic  nozzle  dissemination  (continued) 
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Figure  17:  Aerosii  R-8 12  Dissemination 
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Figure  18;  SipematD-11  Dissemination 
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Several  different  forms  of  titanium  dioxide  were  available  for  testing.  Being  also 
environmentally  benign,  the  Ti02  is  an  attractive  candidate  for  mid-range  size  simulation. 
Primary  particle  sizes  of  the  Tiona  RCL-2  and  R69  (SCM  Chemicals)  and  Tioxide  R-FC2 
(Tioxide,  Inc.)  are  in  the  range  0.2-0.3  micron.  Their  spectra  from  sonic  nozzle  dissemination 
show  median  size  counts  at  0,4-0. S  micron  indicating  reasonably  efficient  de-agglomeration.  It 
should  be  noted  that  these  powders  are  very  hygroscopic  and,  hence,  difficult  to  disperse.  The 
addition  of  vertical  arms  to  the  stirring  rod,  however,  was  sufficient  to  overcome  most  of  these 
problematic  candidates.  The  narrow  spectra  achieved  can  be  useful  especially  as  the  protection 
factor  model  program  indicates  this  size  range  as  particularly  problematic  fhom  the  protection 
aspect. 


Titanium  dioxide  P-2S  is  a  fumed  product  with  an  average  primary  particle  size  of 
21  nm.  Tiosorb  2  (UF02)  is  a  microfine  product  of  Tioxide  with  18  nm  primary  particle  size. 
Unlike  the  hydrophobic  Aerosil  products,  these  are  not  well  de-agglomerated  through  the  sonic 
nozzle  and  are  not  beneficial  candidates  for  further  testing. 

Talc  and  GM  Arizona  Road  Dust  provide  a  useful  broad  spectrum  ranging  from  0.2 
to  7  microns.  Little  difference  was  observed  between  the  spectra  from  the  nozzle  and  that  fhom 
the  fans  only.  The  two  powders  appears  to  exist  originally  as  primary  particles,  whereas  the  nm 
range  primary  particles  of  the  Aerosils  and  Sipemats  experience  more  considerable  adhesive 
forces,  hence  their  agglomerate  natures. 

6.2  Solid  Particle  Settling  Rates 

The  considerable  variances  between  theory  and  measurements  of  solid  aerosol 
settling  rates  noted  in  the  literature  survey  and  the  uncertainty  of  determining  turbulent  energy 
dissipation  parameters  for  stiired  enclosures  necessitated  a  series  of  measurements  specific  to  the 
chamber  itself  Settling  rate  tests  are  carried  out  by  exposing  the  chamber  interior  to  a  continuous 
cloud  of  disseminated  powder  and  then  resealing  it.  This  can  be  done  by  removing  one  perspex 
plate  from  the  window  during  passage  of  the  cloud.  The  mass  concentration  is  monitored  to 
ensure  reasonable  count  rates  before  spectra  are  collected  by  the  LAS-X  each  four  minutes  for  an 
extended  time .  Total  sampling  time  is  dictated  by  the  slowest  settling  rates  of  the  size  bins  of 
interest.  Measurements  can  be  stopped  when  bin  counts  reach  stiitistically  insignificant  levels. 
Various  stirring  scenarios  can  be  examined  by  operating  small  fans  inside  the  chamber.  If  a  true 
quiescent  state,  barring  temperature  differences  within  the  enclosure  space,  is  of  interest,  the 
LAC'X  must  be  situated  outside  the  chamber  and  its  sample  tube  inserted  through  an  opening  in 
the  wall.  It  has  been  observed  that  its  circulation  fan  alone  has  a  significant  effect  on  the  settling 
ratws.  In  all  cases  the  chamber  must  be  completely  sealed  in  order  to  reduce  to  a  minimum 
particle  losses  due  to  air  exchange.  It  should  be  noted  that  results  are  relevant  to  an  enclosure 
with  smooth,  non-conducting  walls.  Effects  such  as  electrostatic  charging  and  rough  surfaces, 
which  may  be  encountered  in  practical  situations,  should  be  considered  for  investigation  as  well. 
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Data  are  analyzed  by  a  semi-logarithmic  optimization  of  the  background-corrected 
particle  counts  in  individual  size  bins  vs.  time  to  yield  best  values  of  the  decay  rates; 

^  =  -3pC(t)  (13) 

In  the  Appendix  Mathcad  utility  programs  needed  for  background  subtraction  of 
LAS-X  data  and  calculation  of  settling  rates  from  this  data  are  described.  In  all  measurements  to 
date  the  decay  rates  in  the  chamber  have  fit  the  single  exponential  function  with  a  correlation 
coefficient  better  than  0.9.  These  results  could  indicate  that  electrostatic  forces  are  not  interfering 
with  the  mechanical  mechanisms  of  settling.  Such  interference  would  be  expected  to  cause  a 
more  rapid  initial  decay  of  small  particles  and  possibly  a  time-dependent  hastening  of  the  larger 
particle  decay.  Since  results  of  settling  measurements  have  been  consistent  over  several  months 
and  under  varying  environmental  conditions,  it  does  not  appear  that  electrostatic  charges  affect 
the  results.  Such  influences  would  be  expected  to  vary  with  changes  in  humidity  at  the  very  least. 

Shown  in  figure  19  are  settling  rates  plotted  against  particle  size  for  talc,  which 
spectrum  provides  a  convenient  tool  for  simultaneous,  polydisperse  measurements.  The  effect  of 
turbulent  mixing  intensity,  increased  between  the  runs  by  operating  a  small  muffin  fkn  in  the 
chamber,  expresses  itself  in  the  settling  rate  curves  which  uniformly  shift  to  higher  levels. 
Measurements  vrith  other,  less  polydispersed  powders  yielded  similar  results  consistent  with 
density  differences.  These  data  can  be  used  directly  in  the  protection  fkctor  model  as  calibration 
parameters.  They  are  useful  also  in  estimating  empirical  values  for  the  coefficient  and  exponent  of 
the  eddy  dissipation  equation  (10),  This  topic  will  be  treated  separately  later. 


6.3  Solid  Aerosol  Protection  Factor  Tests 

Several  series  of  experiments  were  carried  out  with  the  objective  of  investigating  the 
penetration  of  particles  in  the  range  0. 1  -S  microns  into  the  chamber  with  a  ’/entilation  dynamic 
scenario  which  should  force  a  significant  difference  between  the  protection  factors  for  vapor  and 
solid  aerosols.  Their  purpose  is  to  verify  the  results  of  the  model  given  the  same  scenario.  The 
procedures  followed  are  applicable  in  general. 

The  chamber  ventilation  kinetics  are  established  by  first  fixing  the  aperture  width  with 
two  widely  separated  feeler  gauges.  The  gauges  must  remain  in  the  aperture  throughout  the  test 
to  prevent  its  collapse  when  a  negative  pressure  differential  is  applied.  The  perspex  plates  will  be 
sucked  against  the  rubber  seal  and  the  aperture  geometry  may  be  compromised.  For  very  narrow 
apertures  care  must  be  taken  that  the  surftices  are  clean  and  the  plates  are  pressed  together  against 
the  feeler  gauges  before  final  clamping.  The  door  of  the  chamber  must  be  tapes  with  duct  tape 
around  its  entire  perimeter  to  ensure  complete  sealing.  Make  sure  that  all  instruments  are 
operating  and  zeroed  before  doing  this. 

The  desired  negative  pressure  differential  is  established  by  operating  the  blower 
situated  near  the  control  table.  The  mechanical  iris  installed  on  the  air  pipe  must  be  used  to  limit 
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a  -  Spaccromecer  fan  and  muffin  fan  on  in  chamber 
b  -  Spectrometer  fan  on  in  chamber 
c  -  No  fana  in  chamber 
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Figure  19:  Talc  Particle  Settling  Rates 
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suction  on  the  chamber.  Start  with  a  completely  closed  iris  and  with  lowered  voltage  on  the 
Variac  transformer  which  controls  the  blower.  Slowly  raise  the  voltage  back  to  120  volts  and 
open  the  iris  slightly  to  fix  the  desired  negative  pressure  as  indicated  by  the  reading  of  the 
piezoresistive  transducer.  Note  that  for  very  low  pressure  differentials,  the  influence  of  the  'wind' 
impacting  on  the  chamber  is  measurable  and  significant.  Hence,  in  the  case  of  differentials  of 
only  a  few  mmWG  the  final  adjustment  should  be  done  with  the  fans  operating.  Pressure  and 
airflow  should  be  recorded  continuously  during  the  final  adjustment  and  an  average  pressure 
calculated. 


The  air  exchange  rate  can  be  established  either  by  measurement  directly  with  tracer 
gas  (SF5)  or  by  estimation  using  the  measured  airflow  rate  and  assuming  uniform  mixing  of  the 
chamber.  This  assumption  is  generally  valid  when  any  air  mover  is  operating  in  the  chamber  or  a 
temperature  differential  exists  across  the  chamber.  In  borderline  quiescent  tests,  this  assumption 
should  be  verified  through  a  tracer  test. 

Unless  a  specific  exposure  time  scenario  is  dictated  ahead  of  time,  it  may  be 
worthwhile  to  make  a  dry  run  test  with  the  candidate  powder  to  determine  release  and  total 
exposure  time  sufficient  to  produce  statistically  significant  readings  and  proper  instrument  range 
settings.  At  low  air  exchange  rates  and  narrow  aperture  settings  the  penetration  fractions  of  some 
powders  may  be  too  low  for  very  short  exposures. 

Once  all  the  sensors  are  calibrated  and  ready  for  operation,  the  powder  may  be 
released  for  the  chosen  duration,  The  LAS-X  spectrometer,  controlled  by  software  fi'om  the  HP- 
85  computer,  should  begin  accumulating  data  at  the  moment  the  cloud  appears  at  the  chamber,  as 
seen  by  the  RAS-2  sensor  readings.  The  LAS-X  can  record  repeat  spectra  at  a  minimum  4 
minute  interval.  Hence,  if  spectral  data  inside  the  chamber  are  required,  synchronization  with  the 
cloud  will  ensure  capture  of  the  initial  penetrating  particles.  Spectra  should  be  accumulated 
continuously  following  passage  of  the  cloud  and  as  long  as  statistically  significant  data  are 
obtained  from  the  bins  of  interest.  These  data  will  later  provide  size-dependent  particle  decay 
data; 


Rp^R-¥^p 

The  approximate  empirical  protection  factors,  defined  as  the  outside  dose  to  the 
inside  dose  ratio  over  a  given  exposure  time,  are  calculated  ftom  the  analog  mass  concentration 
records.  The  Size-dependent  RAS-2  sensor  responses  are  corrected  by  the  log-normal  parameters 
of  the  challenge  spectrum  and  the  internal  RAM-1  monitor  response  by  the  penetrating  particles 
initial  spectrum  using  the  appropriate  Mathcad  utility  program  (MIE.MCD)  to  be  described  later. 
The  inside  dose  value  is  somewhat  inaccurate  for  long  exposure  times  since  size  composition  of 
measured  mass  concentrations  is  time-dependent,  the  particles  decaying  in  accordance  with  the 
decay  rate  distribution.  Hence,  the  monitor  response  is  time-dependent.  Results  based  on  mass 
concentrations  are  considered  first-order  approximations  unless  the  more  detailed  data  from  the 
spectrometer  are  analyzed  for  their  time-dependent  content. 
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In  one  series  the  same  ventilation  kinetics  were  maintained  while  different  powders 
were  disseminated.  All  disseminations  were  of  0.5  minute  duration,  simulating  the  shortest 
exposure  time  of  interest  in  this  exploratory  series.  The  major  parameters  are  as  follows: 

-  disseminated  mass;  1  gram, 

-  slit  width;  0. 1  mm  (determined  by  feeler  gauges) 

•  slit  depth;  40  mm 

-  pressure  differential;  •!  mmWG 

Results  for  the  first  series  showed  the  expected  phenomena.  Protection  factors 
were  highest  for  the  largest  threat  size  spectra  as  the  large  particles  are  apparently  filtered  through 
the  slit  and  decay  most  rapidly.  These  effects  are  seen  most  clearly  by  comparing  the  plots  of 
Sipemat  Dll  (<0.4  microns)  penetration  with  those  of  talc  (0.3-10  ndcrons)  shotvn  in  figure  20. 
The  upper  pair  of  plots  is  the  Sipemat  test.  The  O.S  minute  exposure  shown  for  two  sensors 
yielded  a  well-developed,  long-lived  contamination  of  the  chamber.  The  Icwer  two  plots  show  a 
series  of  three  consecutive  talc  disseminations,  each  of  O.S  minute,  on  the  same  time  scale  as  the 
plot  above.  The  chamber  concentrations  are  erratic  at  best  and  die  out  very  shortly.  Clearly, 
penetration  of  the  talc,  with  a  mass  median  diameter  about  8  microns.  Is  restricted  in  comparison 
to  the  much  smaller  Sipemat  particles. 

The  follovring  table  gives  the  measured  protection  factors  for  1/2  minute  pulses. 

The  exposure  times  used  in  the  calculations  reflect  the  relative  speeds  of  disappearance  of  the 
m{\jor  mass  fractions  for  each  powder; 


number 

test  powder 

PFp 

exposure  time 
(min.) 

1 

Sipemat  Dll 

16 

20 

*2 

Talc 

70 

1.5 

3 

Arizona  dust 

5 

25 

4 

Talc 

60 

4.5 

5 

A  * 

Ti02  RCL3 

6 

40 

0 

Arizona  dust 

12 

10 

Notes:  *  poor  statistics;  **  10  grams  released  during  10  minutes 

The  higher  protection  factors  measured  for  talc  exposure  could  indicate  both 
considerable  particle  filtering  in  the  aperture  and  rapid  decay  of  the  high  mass  particles.  However, 
the  short  exposures  did  not  permit  the  accumulation  of  sufficient  particle  counts  to  allow 
comparison  of  the  size-dependent  transport  fi'actions  fhom  the  LAS-X  data. 

For  similar  exposure  times,  the  high  protection  factor  for  the  Sipemat  Dll 
compared  to  Arizona  Road  Dust  may  be  a  consequence  of  the  rapid  settling  rate  of  the  smallest 
particles.  Arizona  Road  Dust  has  a  large  portion  of  its  mass  in  the  O.S  to  2  micron  region  in 
which  the  settling  rates  have  been  seen  to  be  the  slowest,  hence  the  protection  factors  are  lowered 
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Figure  20;  Solid  Aerosol  Penetration 
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by  the  ongoing  exposure  to  these  particles,  A  similar  result  occurred  with  the  Ti02  whose 
spectrum  is  in  the  problematic  mid-range. 

The  presentation  of  results  from  the  first  series  is  qualitative  in  nature  since  only  the 
integral  exposure  data  are  being  considered  and  compared  to  expectations  fVom  our  understanding 
of  the  powders  used  in  the  tests.  The  results  will  be  compared  to  model  calculations  in  the  next 
chapter. 


The  second  series  uf  tests  was  designed  to  emphasize  the  size  differentiation  during 
particle  penetration  as  a  function  of  the  differential  pressure  across  the  aperture  using  a  talc 
exposure  of  longer  duration  to  ensure  statistically  sigiiificant  results  .  Three  runs  were  performed 
at  pressure  differentials  of  -1,  >2  and  -3  mmWG  with  the  same  aperture.  Results  of  the 
measurements  are  as  follows: 


Differential  Pressure 

Exposure  Tintes 

Protection 

[mmWG] 

T  [min.] 

t  [min.] 

Factor 

-1 

11 

34 

27 

-2 

11 

29 

15 

-3 

17 

35 

10 

As  expected,  the  protection  factors  measured  fkll  sharply  as  the  pressure  differential 
increases.  This  is  in  part  due  to  the  increased  air  exchange  rate  which  flushes  particles  out  more 
rapidly.  The  extent  to  which  particulate  penetration  itself  was  affected  will  be  investigated 
further.  Note  again  that  at  this  point  only  gross  integral  information  is  provided  by  the  mass 
concentration  sensors.  The  results  will  later  be  compared  to  model  predictions. 

6.4  Solid  Particle  Penetration  Factors 

The  extent  of  particle  filtering  in  the  apertuie  as  a  function  of  pressure  differential 
was  investigated  by  comparing  the  shapes  of  the  particle  size  spectra  of  the  cloud  and  of  the 
penetrating  fraction  for  the  test  at  -  I  mmWG  and  at  -3  mmWG  differential  pressure.  Data  from 
the  LAS-X  spectrometer  provided  counts  of  the  initial  penetrating  fraction.  The  spectrum  of  the 
threat  cloud  was  obtained  by  the  LAS-X  prior  to  the  test  by  sampling  through  the  enclosure  wall. 
In  figure  21  are  the  inside/outside  ratios  of  counts  in  each  size  bin  for  the  two  differential 
pressures.  The  vectors  ratioO.  ratio  1,  and  ratio2  cover  the  range  for  which  differentiation  is 
expected,  between  1  micron  and  the  maximum  size  at  about  7  microns .  The  square  symbols 
show  these  ratio  values  on  the  semi-log  plot.  The  sharp  cutoff  between  0  and  7  microns  at 
-1  mmWG  correlates  well  with  the  model  results  for  the  total  transport  fraction.  Note  that  the 
ratio  at  7. 1  microns  is  statistically  invalid,  being  based  on  single  counts  only.  The  transport 
fractions  at  -3  mmWG  in  the  lower  part  of  Figure  21  show  less  pronounced  filtration,  as  would  be 
expected  from  theoi^.  There  is  no  cutolf  size.  About  40-50%  of  the  pariicles  at  6-7  microns  still 
penetrate  the  aperture. 

The  statistical  fluctuations  of  the  real  data  make  it  difllcult  to  clearly  compare  the 
ratios  for  the  different  tests.  The  solid  linei  shown  are  the  inside/outside  ratios  calculated  from 
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figure  2 1 :  Solid  Particle  Penetration  Factors 
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log  normal  distribution  parameter  estimates  obtained  from  each  data  spectrum.  The  lines  are  not 
extended  further  since  the  theoretical  distributions  are  asymptotic  and  the  real  size  spectra  are,  of 
course,  not.  The  ratios  for  the  test  at  *3  mmWG  pressure  differential  clearly  show  reduced 
filtering  in  comparison  to  the  test  at  -1  mmWO.  These  results  compare  favorably  with  tire  model 
results  shown  in  Figure  1  for  the  same  range  of  pressure  differentials. 

6.5  Aperture  Flow  Characteristics 

Solid  aerosol  dynamics  are  intimately  linked  to  the  ventilation  kinetics  of  the 
enclosure.  The  air  exchange  rate  is  a  major  parameter  controlling  the  particle  balance.  If  the 
assumption  of  ideal  air  mixing  is  valid,  the  equation  R"QA^  can  be  applied,  This,  of  course, 
requires  measurement  of  air  flowrate  and  enclosure  volume.  In  cases  in  which  this  is  not  feasible, 
the  tracer  decay  method,  discussed  earlier,  allows  direct  measurement  of  the  air  exchange  rate. 

The  model  program  is  constructed  to  calculate  the  airflow  directly  ffom  the 
characteristics  of  the  aperture  and  the  pressure  differential  using  equation  5.  This  is  not  essential 
to  the  protection  factor  calculations  if  other  estimates  of  the  air  exchange  rate  are  available.  If  so, 
this  calculation  may  be  bypassed.  It  is  instructive,  however,  to  get  a  feeling  for  flow 
characteristics  by  testing  the  aperture  dynamics  against  theory  and  estimating  the  empirical 
parameters.  This  has  been  done  for  the  40  mm  deep  aperture  at  0. 1  mm  height  with  <1  mmWG 
pressure  through  precise  measurement  of  the  air  flowrate  and  pressure  differentials.  The  flow 
exponent  (n  in  equation  S)  and  the  effective  leakage  area  (and,  hence,  effective  height)  were  then 
estimated  by  a  Mathcad  utility  program  designed  for  this  purpose.  A  sample  run  of  the  program 
for  this  calculation  is  given  in  the  Appendix. 

Pressure  differentials  in  the  range  0.4  to  0.8  mmWG  were  measured  with  the 
Endevco  piezoresistive  transducer  and  recorded  continuously  on  the  strip-chart  recorder.  Air 
flowrate  at  the  blower  output  was  measured  with  the  Gilian  Instrument  Corp.  Gilibrator  Primary 
Flow  Calibrator.  The  Gilibrator  can  handle  flows  up  to  30 1pm  only,  hence  limiting  the  range  of 
application.  However,  its  being  a  primary  calibration  gave  the  advantage  of  reduced  uncertainty. 

Estimation  of  the  best -fit  parameters  for  Equation  5  using  five  data  points  yielded  a 
correlation  coefficient  of  0.993.  The  flow  exponent  value  for  the  nominal  0. 1  mm  aperture  is  0.52 
and  the  effective  height  of  the  aperture  is  0.15  nun.  The  difference  between  the  nominal  and  the 
effective  aperture  height  is  reasonable,  considering  the  roughness  of  the  height  measurement  using 
only  feeler  gauges.  Extrapolation  of  the  characteristic  curve  to  -1  mmWG  yielded  a  flowrate  of 
29  1pm  which  was  also  observed  by  direct  measurement 
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7.  SOLE)  AEROSOL  PROTECTION  FACTOR  MODEL  CALIBRATION  AND 

VERIFICATION 

7.1  Model  Calibration 

Results  of  the  model  program  given  earlier  were  based  solely  on  theoretical 
calculations  of  the  parameters  determining  solid  aerosol  behavior  in  the  enclosure.  These  include 
the  air  exchange  rate  R,  the  size-dependent  settling  rate  /S^and  filtration  factor  fp .  Calibration 

of  the  model  prior  to  verifying  it  against  solid  aerosol  protection  fhctor  measurements  requires 
consideration  of  the  validity  of  calculated  parameter  values. 

As  stated,  the  air  exchange  rate  can  either  be  calculated  theoretically  or  estimated 
from  measurement.  The  current  version  of  the  model  uses  the  crack  flow  equation  and,  hence, 
empirical  values  of  the  relevant  parameters  (effective  height  and  flow  characteristic  exponent,  n) 
have  been  used  to  calibrate  this  segment  of  the  model  since  their  validity  has  been  assured. 

The  theoretical  model  for  the  size-dependent  particle  settling  rates  (Equation  8)  does 
not  produce  a  satisfactory  fit  to  the  experimental  data  for  the  enclosure  (the  rate  of  change  of 
settling  rates  with  particle  size  for  particles  larger  than  1  micron  is  lower  experimentally  than  the 
theoretical  value)  although  it  does  reproduce  published  results  of  other  researchers  who  have  used 
it  to  fit  their  empirical  data.  Attempts  to  find  best-fit  parameter  values  using  the  Asystant 
scientific  programming  language's  curve-fitting  algorithm  found  the  problem  to  ha  rather  ill- 
conditioned.  The  settling  rate  model  was  programmed  as  a  nested  user  function  and  run  through 
the  optimization  algorithm  with  many  initial  parameter  values.  The  solutions,  graphs  of  which  are 
shown  in  Appendix  C.4  along  with  the  functions,  showed  the  fitting  function  to  be  too  inrensitive 
to  the  turbulent  energy  intensity  to  be  of  value.  A  sampling  of  the  results  is  given  here; 

Turbulent  intensity  (ke)  Exponent  (n)  Goodness  of  fit  (R^) 

1000  2,038  0,9335 

5000  2.222  0,9359 

10000  2.301  0,9361 

Since  no  reasonable  combination  of  parameters  succeeded  in  reproducing  the 
observed  settling  rates  at  higher  sizes,  empirical  data  were  adopted  to  calibrate  the  model.  The 
form  of  the  size-dependency  was  found  to  be  identical  for  several  powders  tested  to  date  so  only  a 
weighting  fkctor  was  seen  to  be  required  in  the  model  to  adjust  the  calculation  for  diflerent 
powder  densities. 

The  transport  of  particles  through  the  aperture  appears  to  be  modeled  well  enough 
at  this  point  to  permit  use  of  the  theoretical  equations  in  the  model.  Of  course,  fUrther  runs  with 
additional  powder  challenges  are  recommended,  but  there  is  no  objective  reason  to  reject  the 
model  as  it  stands  in  view  of  the  positive  results  of  the  penetration  tests  at  low  differential 
pressures. 
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7,2  Model  Output  Verification 

The  model  was  run  with  the  calibrations  described  in  the  previous  section  included 
and  with  parameters  identical  to  those  of  the  test  experiments  described  previously,  including  the 
actual  exposure  times  used  to  calculate  the  protection  factors  from  the  mass  concentrations.  The 
model  results  for  both  particles  and  vapors  are  compared  to  the  experimental  results  in  the 
following  table.  Recall  that  these  runs  were  designed  to  accentuate  the  difference  between  the 
two  protection  factors  by  using  minimum  pressure  differential  (-1  mmWG)  and  a  very  narrow  (0. 1 
mm)  and  deep  (40  mm)  aperture.  The  model  calculations  in  the  table  include  log-normal 
weighting  with  the  known  challenge  spectrum  in  each  case  and  the  mass  concentration  ratios 
comprising  the  measured  PFp  include  sensor  corrections  in  accordance  with  the  same  spectra.. 

Model  Results 


# 

Test  Powder 

‘*g 

[microns] 

Exposure  Time 
[mins.] 

PFp(measured) 

PFp(solid) 

PFv(vapor) 

1 

Sipemat 

0.18 

20 

16 

13 

8 

*2 

Talc 

1,7 

1.5 

70 

140 

114 

3 

Arizona  dust 

0.7 

25 

5 

10 

6 

4 

Talc 

1.7 

4.5 

60 

50 

34 

5 

Titanium  dioxide  0.6 

40 

6 

6 

4 

*•6 

Arizona  dust 

0.7 

10 

12 

21 

16 

Notes:  *  poor  statistics;  **  10  grams  during  10  minutes 

The  numerical  comparison  between  measured  and  calculated  values  of  PFp  is 
reasonably  good  considering  the  brevity  of  the  experiments.  The  0,5  minute  pulses  led  to  poorer 
statistics  than  could  be  obtained  with  more  extended  exposures.  However  they  do  show  that  the 
experimental  system  can  handle  the  desired  span. 

Note  that  the  differences  between  the  vapor  and  particulate  protection  factors  are 
also  as  expected  by  the  model.  The  smallest  difference  is  for  titanium  dioxide,  which  has  a  narrow 
size  spectrum  in  the  0.3-1  micron  range.  From  Figure  3  we  would  expect  the  least  difference 
between  the  factors  in  this  region.  Sipemat,  with  an  narrow  spectrum  near  0.2  microns,  produced 
a  factor  2  difference  between  the  solid  and  the  vapor  protection  factors,  again  as  expected  fhom 
the  model  results.  Also  according  to  expectations,  the  broad  Arizona  road  dust  spectrum  below  1 
micron  did  not  produce  a  significant  difference.  A  most  clear  difference  occurred  with  the  talc 
exposure  in  test  4.  The  talc  spectrum  includes  a  significant  portion  of  particles  over  1  micron  and, 
hence,  according  to  the  model  results,  a  large  difference  between  PFp  and  PFy  is  expected.  The 
talc  run  in  test  2  does  not  appear  to  be  representative  as  the  measured  exposure  inside  the 
enclosure  was  too  short  to  ensure  statistically  significant  results.  For  this  reason  the  next  set  of 
runs  was  done  with  longer  challenge  duration  to  ensure  a  fUlly  developed  inside  exposure.  Note 
that  these  brief  experiments  were  carried  out  far  from  the  equilibrium  exposure  levels.  At  the  low 
exchange  rates  used,  equilibrium  exposure  levels  would  have  been  reached  only  after  several 
hours. 
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The  next  set  of  experiments  investigated  the  influence  of  the  pressure  differential  in 
determining  the  particle  protection  factor,  Three  runs  were  carried  out  with  talc,  the  challenge 
with  the  largest  expected  differences  in  the  factors.  Three  pressure  differentials,  again  in  the  very 
low  rarvge  expected  to  still  show  some  size  differentiation,  were  used.  Challenge  durations  were 
extended  sufficiently  to  effect  a  ililly-developed  inside  exposure.  The  previously  shown  results  of 
these  tests,  along  with  the  model  calculations  are  given  in  the  following  table. 


Model  Results 

Pressure 

Exposure  Time 

PFp(measurcd) 

PFp(solid) 

PFv(vapor) 

[mmWG] 

[min,] 

-] 

11 

27 

IS 

5.5 

-2 

11 

15 

15 

7 

-3 

17 

10 

12 

3.6 

The  model  results  compare  reasonably  well  with  experiment  for  the  solid  aerosol 
protection  factors,  keeping  in  mind  the  uncertainty  due  to  the  time  dependency  of  the  sensor 
responses.  The  decrease  in  vapor  protection  with  increasing  pressure  differendal  is  obvious  as  the 
air  exchange  rate  increases.  The  solid  aerosol  protection  factors  tracked  the  expected  decrease  as 
thu  increasing  pressure  differential  pulls  more  and  larger  particles  into  the  enclosure. 

7,3  Sensitivity  Analysis 

Having  seen  that  a  reasonable  correlation  between  the  expected  and  the  measured 
results  is  obtained,  the  model  may  be  then  used  for  parameter  sensitivity  testing.  Such  numerical 
testing  helps,  for  example,  to  determine  under  what  circumstances  a  significant  difference  between 
the  vapor  and  the  solid  aerosol  protection  factors  may  occur  or  what  must  be  done  to  prevent 
particles  with  a  given  size  distribution  from  penetrating  an  enclosure. 

In  the  model  runs  used  for  Figure  22,  all  the  parameters  are  those  previously  used  in 
the  verification  runs  for  a  talc  challenge.  Only  the  aperture  height  is  changed  between  the 
diagrams.  In  each  diagram  three  sets  of  protection  factors  are  plotted  against  the  differential 
pressure:  vapor,  solid  aerosol  with  quiescent  mixing  and  solid  aerosol  with  more  turbulent  mixing 
conditions  (two  small  fans  operating  in  the  chamber).  As  seen  in  the  experimental  results,  a  O.IS 
mm  slit  width  with  turbulent  mixing  particularly  improves  the  solid  aerosol  protection  at  all 
pressure  differentials  as  large  particles  are  filtered  out  during  penetration  and  the  remainder  decay 
to  the  surfaces  more  rapidly  .  This  influence  of  turbulent  impaction  is  sharply  accentuated  in  the 
curve  for  quiescent  mixing  as  the  solid  aerosol  protection  factors  approach  those  of  the  vapor. 

The  decay  rates  of  the  particles  have  slowed  drastically  and  they  are  flushed  out  more  nearly  at 
the  rate  of  the  vapor  molecules. 

With  a  1.0  mm  slit  width  the  filtering  of  particles  is  already  insignificant  and  the  solid 
aerosol  protection  factors  are  rapidly  approaching  those  of  the  vapor,  being  kept  higher  primarily 
by  turbulent  impaction.  With  a  2.0  mm  slit  width  the  difference  between  the  vapor  and  the 
particulate  protection  factors  is  inconsequential  for  all  practical  purposes. 
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8,  CONCLUSIONS 

The  primary  goal  of  the  research  program  was  to  develop  the  methodology  needed 
to  determine  solid  aerosol  protection  factors  for  enclosures  and  to  express  the  controlling 
mechanisms  in  model  calculations. 

A  model  program  has  been  constructed  which  calculates  theoretical  vapor  and  size> 
dependent  solid  aerosol  protection  factors  for  an  enclosure  with  pressure-driven  penetration.  The 
model  integrates  the  most  current  studies  of  particulate  penetration  factors  and  settling  rates  into 
the  theoretical  equation  for  the  protection  factor  of  a  well-mixed  enclosure.  The  model  program 
is  \vritten  in  the  Mathcad  scientific  language,  providing  an  interactive  tool  appropriate  for 
protection  estimation,  parameter  sensitivity  studies,  and  test  procedure  design. 

An  experimental  system  has  been  assembled  which  can  provide  real  data  on  the 
parameters  relevant  to  the  calculation  of  solid  aerosol  protection  factors  for  enclosures  The 
system  can  be  used  to  produce  data  for  model  calibration,  verification  and  prediction  as  well  as 
for  constmction  of  a  database  of  generally  applicable  parameter  values.  Together  with  the  model 
program  the  facility  will  aid  in  optimizing  integrity  testing  options  for  shelters  and  vehicles. 

When  combined  with  empirical  data,  model  results  to  date  correlate  well  with 
observed  behavior.  It  appears  that  size-dependent  filtration  of  particles  in  apertures  will  be 
limited  to  the  millimeter  width  level  at  pressure  differentials  up  to  only  several  millimeters  water 
gauge.  The  extent  of  turbulent  impaction  of  particles  within  the  enclosure  will  strongly  affect  the 
particulate  protection  factors.  However,  fUrther  work  is  needed  to  provide  better  characterization 
of  particle  settling  in  a  stirred  chamber  in  terms  of  the  turbulent  energy  content  before  the  effect 
can  be  satisfactorily  quantified.  It  must  be  stressed  that  the  sensitivity  of  the  model  output  to  the 
settling  rates  and  the  wide  variance  in  the  turbulent  dissipation  parameter  values  seen  in  the 
literature  warrants  caution  in  their  application.  Filtration  of  specific  particles  in  rough  apertures, 
phenomena  such  as  possible  electrostatic  charges,  surface  effects,  resuspension,  and 
agglomeration  will  have  to  be  characteri.zed  if  the  phenomena  determining  the  protection  factors 
for  various  challenges  are  to  be  fully  understood.. 


9.  RECOMMENDATIONS 

The  limited  timeframe  of  the  program  and  the  onset  of  inclement  weather  precluded 
implementation  of  a  comprehensive  series  of  tests.  Therefore,  the  tests  performed  were  designed 
to  accentuate  the  mechanisms  which  should  lead  to  differentiation  between  the  vapor  and  the 
particulate  protection  factors.  This  having  been  done  and  the  results  having  demonstrated  the 
ability  of  the  system  to  provide  quantitatively  significant  infonnation,  it  is  felt  that  there  is  now  a 
basis  for  specific  threat-oriented  questions  to  be  treated. 

The  analytical  scheme  adopted  for  the  exploratory  phase  of  the  program  has  made 
use  of  direct  particle  size  count  dnta  and  integrated  mass  concentration  data.  It  does  not  address 
other  ways  of  looking  at  particulate  data,  such  as  the  mass  spectral  distribution  of  polydisperse 
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threats.  It  is  assumed  that,  further  along,  those  concerned  with  specific  threat  scenarios  will  adopt 
the  most  appropriate  definitions  and  make  the  conceptual  link  to  the  particulate  protection  factor. 

The  model  is  based  on  the  assumption  that  the  settling  process  is  a  perfect  sink. 
Particles  reaching  the  walls  are  lost  to  the  system  and  do  not  pose  a  further  threat.  Thought 
should  be  given  to  more  practical  aspects  of  particles  in  a  confined  enclosure,  such  as  in  a  crew 
compartment,  which  may  be  mechanically  resuspended  and  pose  a  recurrent  threat. 

The  shift  from  studying  a  well-defined  enclosure/aperture  system  to  testing  an  entire 
vehicle  or  shelter  will  have  to  be  done  with  caution.  So  far,  candidate  powders  of  opportunity 
were  used  to  help  accentuate  size-dependent  phenomena.  Calibrated  monodisperse  particles 
should  be  adopted  also  to  reduce  the  complexity  of  the  mass  concentration  sensor  responses, 
especially  in  situations  where  experimental  uncertainty  is  already  high  due  to  the  gross  integral 
nature  of  a  whole-vehicle  test.  Considerable  care  will  have  to  be  taken  also  to  ensure  uniform 
exposure  of  a  vehicle  to  the  threat  cloud  in  static  scenarios,  preventing  erroneous  conclusions  due 
to  spatial  variances  in  the  test  concentration.  Even  more  critical  will  be  the  need  to  produce  and 
characterize  realistic  exposure  profiles  in  complex  dynamic  scenarios,  such  as  moving  vehicles 
traversing  a  transient  particulate  cloud. 

In  order  to  provide  a  basic  heuristic  tool,  the  model  constructed  for  this  program 
involved  the  most  simplistic  mass  balance  assumptions  and  the  experimental  setup  upon  which  its 
verification  is  based  simulates  this.  It  cannot  be  emphasized  too  strottgly  that  the  solid  aerosol 
dynamics  are  driven  first  by  the  ventilation  kinetics  of  the  system  These  should  always  be 
delineated  by  means  of  ei..pirical  tracer  gas  methods  before  the  air/particle/enclosure  system  can 
be  considered  ftilly  characterized. 

It  was  the  intention  at  the  outset  of  this  program  to  assemble  a  broad,  flexible,  state 
of  the  art  instrumentation  and  computational  system  capable  of  independently  providing  reliable 
data  and  analyses  to  answer  as  many  questions  as  may  arise  in  the  future  regarding  the  protection 
of  vehicles  and  enclosures  from  both  vapor  and  solid  aerosol  challenges.  It  is  the  intention  of  this 
Technical  Report  and  its  ^  ppendices  to  aid  the  next  users  of  the  system  in  its  successful 
application.  The  key  words  here  are  likely  to  be  flexibility  and  interaction.  Instnamental 
components  should  be  configured  to  answer  specific  questions  and  reconfigured  again  when 
results  raise  further  questions,  as  so  often  happens.  Model  and  utility  programs  developed  so  far 
should  be  altered,  upgraded  or  discarded  in  accordance  with  their  usefulness.  They  are  designed 
to  be  changed  readily  and  deleted  just  as  readily.  However,  the  user  Is  cautioned  not  to  rely  only 
on  what  is  offered  here,  rather  to  familiarize  himself  thoroughly  with  the  instrument  manuals, 
programming  codes  and  literature  reference;,  as  well.  Analysis  of  solid  aerosol  data  can  often  be 
complex  at  least,  open  to  subjective  interpretation  at  best,  and  downright  tricky  at  worst.  Hence, 
the  importance  in  applying  the  right  methods  to  the  right  questions.  It  is  hoped  that  some 
groundwork  for  flexible  responses  has  been  achieved  here. 
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APPENDIX  A 

Protection  Factors  for  Vapor  and  Solid  Aerosols 


The  following  provides  a  brief  and  simplified  derivation  of  equations  (1)  and  (2)  in 
the  text  which  define  the  protection  fkctors  for  vapor  and  solid  aerosol  penetration  into  leaky 
enclosures.  The  major  addition  here  is  that  the  exposure  resulting  from  continued  occupation  of 
the  enclosure  after  the  outside  challenge  has  passed  is  taken  into  account. 

The  analysis  assumes  spatially  uniform  exposure  of  all  leakage  areas  to  the 
challenge  concentration  outside.  It  is  further  assumed  that  the  infiltrating  air  carries  with  it  the 
true  outside  concentration  of  ciiallenge  material  at  the  enclosure  boundary.  This  ideally  ignores 
any  possible  boundary  layer  effects  which  could  arise  under  various  circumstances.  The  vapor  is 
non>condensing  and  non-depositing. 

The  enclosure  is  characterized  by  an  infiltration  rate,  R,  in  units  of  air 
changes/hour,  determined  by  the  ventilation  dynamic  scenario  (wind  speed  and  stability,  thermal 
gradients,  leakage  areas,  etc).  Assume  the  enclosure  to  be  exposed  to  a  "square"  cloud  of 
contaminant  in  which  the  concentration,  CE,  increases  abruptly  to  a  value  CO  and  then  decreases 
abruptly  to  zero  after  flrtite  time  T.  The  rate  of  change  of  the  inside  concentration  will  be; 

This  has  the  solution,  for  0  <  t  <  T  when  CE“CO: 

C/i(0  =  CO'(l-exp(-.R.7’)) 

and,  for  t  >  T  when  CE"0; 

aj(<)=.c/{r)-exp(-R.(r-r)) 

Hence: 

C/j (/)  =  CO  •  (exp(i? ■  r)  - 1)  ■  exp(-/?  l) 

We  require  the  concentration-time  integrals  in  order  to  calculate  the  dosages  outside 
and  inside  the  enclosure.  At  time  t-T  the  outside  dose  is;  DO  =  CO  - T 


Inside  the  enclosure,  for  0  <  t  <  T; 


D/i  =  jc/j(r)t// 

0 

J3/,=^.(R.7--(l-«p(-R.r))) 


For  t  >  T: 
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1 

Dl^  =  \CI^{t)dt 
T 

Hence; 

exp(7?  •  r))  •  (exp(-7?  •  /)  -  exp(-i?  •  /)) 

A 

The  protection  factor  is  defined  as  the  ratio  of  the  dose  which  would  result  A-oni 
exposure  to  the  outside  concentration  to  the  dose  accumulating  in  the  enclosure  up  to  time  t  >  T; 

PF» - 

D/j+JD/j 

Hence; 

pp _ ^ _ ^ 

V  “  exp(-/?  •  r)) + ( 1  -  exp( « •  f))  •  (exp(-/?  •  /)  -  exp(-;?  •  f)) 

Reducing  this  yields  equation  (1); 

PP _ R'T 

V  "  .«.r+exp{-/?-/).(l-exp(/?T)) 

Note  the  special  case  in  which  the  cloud  passage  time  equals  exactly  the  inside 
exposure  time,  i.e.,  the  occupants  exit  the  enclosure  immediately; 

/?*  T 

""  /i.r-fi-expi-z^-r)} 

Addressing  now  a  solid  aerosol  challenge,  two  effects  that  influence  the  interior 
dosage  and,  hence,  the  derivation  of  the  protection  factor,  are  the  filtering  of  particles  as  they 
pass  through  the  enclosure's  leakage  areas  (cracks,  seams,  apertures,  outlets,  etc.)  and  the 
deposition  of  particles  on  interior  surface  areas.  Resuspension  of  particles  can  also  conceivably 
affect  the  overall  balance  at  a  later  time  as  can  evaporation  of  condensed  vapors,  but  is  not 
considered  here.  Both  phenomena  are  particle-size  dependent  and  as  such  cannot  be 
characterized  completely  for  a  challenge  scenario  on  the  basis  of  ventilation  dynamics  alone 
without  knowing  the  characteristics  of  the  challenge.  We  look  at  the  basic  balance  equation  for 
the  enclosure  with  the  filtering  factor,  fp,  which  reduces  the  penetrating  particle  concentration, 
and  the  decay  rate,  /3^,  which  contributes  to  the  removal  of  particles; 

= 4  ■R-[co,(i)  -  a,{i))  -  /J,  •  c/,(() 

Wc  can  combine  the  two  removal  terms  into  one  effective  term  (see  Engelmann, 

1992a); 
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Hence: 


^C/,(<)  =  /,  R-CO^O) - R,  ■  Cl,{l) 


Solving  the  equation  of  the  form 


yields; 


a  +  b'  X 


?(-R,-l) 


with  the  given  boundaiy  conditions  for  0<t<T: 


CV„W  =  ^-CO,.(l-.xp(-«,.2)) 


Solving  for  t>T: 


a.(')  =  =^'CO,-(l-«ip(-/e,'r))’Cxp(-J!,.(/-7-)) 


Again,  we  require  doses  for  the  two  time  periods; 
D/„=|c/,,(rVr  and  i 


fn-R  (  1  exp(-/?,‘r)'i 


O/,.  =  ^ •  CO,  •  (l - «p(ii, .  7-))  i-. («p(-«, .()- exp(-i!, .  r)) 


By  definition  we  have; 

p2 

Combining  and  reducing  terms  yields  equation  (2); 


where  DO^^CO,>T 


PF  = 

f 


R.‘T 


^^{r,  •  7’+  exp(-;?, . /)-(l -  exp(/?,  •  7'))) 
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APPENDIX  B 


Protection  Factor  Modei  Program 

The  model  program  for  protection  factor  calculation  is  written  in  the  Mathcad 
programming  language.  The  language  combines  equations,  graphics  and  text  in  the  same 
document.  Equations  are  written  directly  into  the  program  as  they  would  be  written  by  hand  and 
are  solved  sequentially.  Text  and  graphs  are  included  and  results  are  displayed  as  desired  so  that 
the  actual  program  has  the  form  of  a  report.  The  model  program  is  included  liet  e  in  the  following 
pages.  It  is  assumed  that  the  reader  has  familiarized  himself  with  the  language.  Only  the 
minimum  amount  of  text  necessary  to  follow  the  calculations  is  included  within  the  body  of  the 
model.  In  order  to  preserve  its  form,  additional  explanations  and  comments  are  given  here, 
followed  by  the  model  itself 

Pagel:  PARAMETER  VALUES 

All  the  parameters  relevant  to  the  model  run  are  concentrated  on  the  first  page. 
These  include  those  needed  for  particle  transport  (density,  viscosity,  temperature  and  physical 
constants),  modelling  constants  (particle  size  range,  indices  and  unit  multipliers),  system 
parameters  (enclosure  dimensions,  aperture  dimensions,  differential  pressures  and  flow  constants) 
and  experimental  parameters  (turbulent  dissipation  parameters,  empirical  settling  rates,  cloud 
duration,  exposure  time,  challenge  spectral  parameters  and  an  optional  challenge  spectrum).  The 
'current  test  ID*  is  revised  to  indicate  the  test  fhom  which  the  data  were  taken.  Of  course,  all 
parameters  listed  on  the  page  can  be  altered  for  any  run  as  long  as  their  values  reflect  physical 
reality  (recall  the  GIGO  Syndrome), 

Page  2:  PRELIMINARY  CALCULATIONS 

Air  flowrate,  required  for  calculation  of  vapor  and  particulate  transport,  is  calculated 
fhom  the  'crack*  equation  with  the  differential  pressure  as  the  driving  force.  Note  that  differential 
pressure  is  an  independent  variable  given  five  values  from  1  to  3  mmW(]i  in  this  run.  Hence  all 
calculations  depending  on  it  are  vectors  with  the  index).  All  particle  size-dependent  calculations 
are  indexed  with  i,  the  range  of  sizes  from  0.1  to  10  microns  defined  at  the  beginning  of  the 
program.  Cunningham's  correction,  terminal  velocities,  diffusion  coefficients,  settling  parameters 
and  the  lognormal  distribution  function  are  calculated  here  for  all  future  use.  The  Reynold's 
numbers  for  flow  through  the  aperture  are  printed  for  reference,  although  not  included  in  any 
calculations. 

PENETRATION  TRANSPORT  -  DEPOSITION  BY  MOLECULAR  DIFFUSION 

The  empirical  equations  for  transport  through  a  rectangular  channel  mentioned  in  the 
literature  search  (Schwendiman  and  Sutter,  1977)  are  adopted  for  use  heie.  In  order  to  simplify 
the  programming,  if  the  condition  /3  -<  0. 1  exists,  and  the  alternutive  expression  given  should  be 
used,  the  change  is  done  manually  later  on.  Values  of  this  parameter  along  with  selected  values  of 
the  difRision  coefficient  and  the  partial  transport  fraction  are  given  on  page  3  against  paiticle  size 
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mQ.  These  results  should  be  examined  at  the  first  run  with  new  parameters  to  determine  if  the 
alternative  expression  should  be  used. 

Page  3:  PENETRATION  TRANSPORT  -  TLTRBULENT  IMPACTION 

As  stated  in  the  report,  deposition  velocities  for  narrow  apertures  and  consideration 
of  rough  surfaces  have  not  been  treated  here.  The  terminal  velocities  are  used  for  now,  resulting 
in  minimum  wall  losses.  Selected  values  of  the  results  are  printed  for  reference. 

Page  4;  PENETRATION  TRANSPORT  -  GRAVITATIONAL  SETTLING 

The  equation  given  in  Schwendiman  and  Sutter  (1977)  is  shown  along  with  selected 

results. 


TOTAL  TRANSPORT  FRACTION 

The  three  mechanisms  are  combined  into  a  final  expression  for  the  fVaction  of 
particles  entering  which  are  transported  through  the  aperture.  Note  that  the  small  square  symbol 
afier  an  equation  written  in  Mathcad  indicates  that  the  expression  is  not  included  in  the 
calculations,  rather  it  is  part  of  the  text.  The  equation  is  essentially  'turned  off.  At  this  point  the 
user  must  chose  the  proper  expression  dependent  on  the  value  of  ^ ,  The  correct  expression  is 
'turned  on'  and  the  alternate  is  'turned  off. 

Page  5;  TOTAL  PENETRATION  TRANSPORT  MATRIX 

1  he  total  penetrating  fractions  are  shown  in  a  two-dimensional  matrix  as  a  function  of 
particle  sizes  (rows)  and  differential  pressures  (columns).  The  particle  size  range  is  shown  at  the 
lefi  for  reference.  Pressure  values  increase  to  the  right  in  the  matrix. 

As  noted  in  the  comment,  negative  values  of  the  transport  traction  at  the  highest  sizes 
and  lowest  pressures  are  liable  to  appear  due  to  exaggerated  losses  by  the  combined  mechanisms. 
This  is,  of  course,  not  physically  realized  but  rather  due  to  the  simplicity  of  the  model.  More 
detailed  modelling  of  the  deposition  velocities  in  narrow  crack  flow  are  needed  but  is  beyond  the 
current  scope  of  the  program.  The  negative  valued  results  are  removed  later. 

The  plot  shows  the  transport  fractions  against  particle  size  over  the  entire  size  range 
with  the  differential  pressure  as  the  parameter  ( pressure  increases  with  the  higher  curves) . 

Page  6:  STIRRED  SETTLING  OF  PARTICLES 

The  particle  settling  rates  are  calculated  here  using  ther  equation  given  by  Chen,  et  al 
(1992).  The  two  parameters,  ke  and  ne,  define  the  turbulent  dissipation  intensity.  The  resulting 
vector  of  theoretical  settling  rates  (solid  line)  is  plotted  against  particle  size  along  with  the 
empirical  values  supplied  initially  (dotted  line).  The  extent  of  the  goodness  of  fit  for  the 
theoretical  values  is  clearly  seen.  This  submodel  may  be  used  to  manually  optimize  the  two 
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p.irameters  or  to  analyze  their  sensitivity,  Automatic  optimization  was  done  separately  with  the 
ASYSTANT  language.  The  best  fit  parameters  have  been  used  here.  As  is  readily  seen  in  the 
plot,  the  fit  is  not  satisfactory  in  this  instance.  The  size  of  minimum  settling  rate  seen 
experimentally  (O.SS  microns)  is  reproduced  and  the  settling  rate  gradient  with  increasing  size  is 
too  great.  No  combination  of  parameter  values  was  able  to  reproduce  these  results  within  the 
framework  of  the  chosen  model. 

For  this  reason,  empirical  values  are  used  in  parallel  with  theoretical  values  for 
calculating  the  protection  factors  for  solid  aerosols.  At  some  time  it  will  have  to  be  decided 
whether  the  sensitivity  to  the  differences  in  the  settling  rates  warrants  tiirther  theoretical  work  on 
improving  the  published  model  for  the  scenarios  to  be  researched  at  a  later  date. 

Primarily  as  a  tool  for  comparison  of  integral  mass  concentration  data  with  theoretical 
results,  settling  rates  have  been  weighted  with  the  lognoitnal  distribution  of  the  challenge  size 
spectrum  The  proximity  of  the  resulting  weighted  numbers  shows  that  in  this  case  the 
discrepancies  may  not  be  too  severe.  This  aspect  should  be  given  attention  in  all  cases  where 
uncertainty  is  high  before  using  the  theoretical  results  further. 

Page  7:  PROTECTION  FACTORS  FOR  THE  EXPOSURE  SCENARIO 

The  vapor  protection  factor  is  based  on  the  air  exchange  rate,  Rj,  for  each  pressure 
diiferential  value.  This  is  calculated  here  based  on  the  assumption  of  complete  mixing  in  the 
enclosure.  This  assumptiom  has  been  verified  for  the  enclosure  under  like  ventilation  kinetics  by 
the  use  of  tracer  gas.  The  flow  rate  calculated  by  the  crack  flow  equation  and  displ?y'''<  here  in 
units  of  m^/hr  and  1/min  has  been  calibrated  by  direct  measurement  with  a  primary  ?■.,  ■  ■  tandard. 

The  total  particulate  removal  rate,  Rpy,  includes  the  air  exchange  rate  and  the 
settling  rate.  Representative  values  at  0.1,  1.0  and  10  microns  are  displayed  and  can  be  compared 
to  the  vapor  removal  rate. 

The  protection  factors  are  calculated  next  using  the  equations  derived  in  Appendix  1 . 
Negative  results  for  the  particulate  protection  factor  due  to  exaggerated  filtration  of  large 
particles  at  low  pressure  differentials  are  eliminated  by  setting  them  equal  to  the  last  positive 
result,  making  this  occurence  obvious  in  the  final  result  tables. 

Page  8:  THEORETICAL  PROTECTION  FACTORS 

This  page  shows  a  summary  of  the  relevant  parameters  for  the  scenario  and  tables  of 
the  results.  The  vapor  protection  factors,  Pv,  are  shown  as  a  vector  of  results  in  order  of 
ascending  pressure  values.  The  particulate  protection  factors,  Pp,  are  shown  as  a  two- 
dimensional  matrix  in  particle  sizes  (rows)  and  presusres  (columns).  All  the  factors  are  displayed 
on  the  plot  below.  The  vapor  factors  are  straight  lines,  being  independent  of  particle  size.  The 
influence  of  the  particle-sp'icific  mechanisms  are  seen  clearly  at  the  lowest  and  highest  size  ranges. 
The  y-axis  of  the  plot  has  been  cut  off  at  50  in  order  emphasize  the  efiects  at  the  lower  values. 
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Page  9:  WEIGHTED  PROTECTION  FACTORS  FOR  THE  CHOSEN  SIZE  SPECTRUM 

On  the  first  page  of  the  program,  the  option  was  given  to  provide  a  particle  size 
spectrum,  Nj,  of  any  desired  spectral  distribution.  This  spectrum  is  now  used  to  weigh  the 
theoretical  size-dependent  protection  factors  accordingly.  This  weighting  yields  a  single 
numerical  value  for  the  given  challenge.  These  values,  wPFj,  are  plotted  against  the  pressure 
differentials  together  with  the  corresponding  vapor  protection  factors. 

Page  10:  WEIGHTED  PROTECTION  FACTORS  FOR  LOGNORMAL  DISTRIBUTION 

If  the  challenge  size  spectrum  has  been  characterized  by  a  lognormal  distribution, 
whose  parameters  have  been  supplied  at  the  beginning  of  the  program,  the  results  can  be  weighted 
by  this  optimum  distribution.  The  distribution  is  shown  in  the  small  plot,  A  table  of  optional 
statistical  definitions  (Reist,  1984)  is  provided  as  a  tool  for  visualizing  other  characteristics  which 
may  be  of  specific  interest.  Any  definition,  dp,  can  be  calculated  by  chosing  the  corresponding 
parameter,  p,  from  the  table. 

Since  the  theoretical  distribution  is  asymptotic  and  the  real  spectrum  is  not,  the 
weighting  Rmction  is  truncated  to  the  end  of  the  real  spectrum  by  cutting  off  the  index  i.  The 
protection  factors  are  then  weighted  by  the  truncated  distribution  and  diplayed  in  the  plot  along 
with  the  vapor  protection  factors, 

Pages  11-12  ,  PROTECTION  FACTORS  USING  EMPIRICAL  SETTLING  DATA 

Since  the  theoretical  settling  rate  results  may  not  be  satisfactory,  the  protection 
factors  arc  recalculated  here  using  the  empirical  data  supplied  at  the  beginning.  The  data  are 
plotted  again  for  reference.  The  weighting  factor,  mult,  is  used  to  chose  the  correct  level  of 
turbulent  dissipation  for  which  the  data  are  representative  ( see  figure  18).  The  protection  factors 
are  again  weighted  by  the  lognormal  distribution  and  displayed  on  page  12  with  the  theoretical 
factors  and  the  vapor  factors.  Here  >he  effect  of  the  discrepancies  between  the  empirical  and  the 
modelled  settling  ra^es  is  clearly  seen.  The  decision  can  be  taken  at  this  point  as  to  whether  more 
effort  should  be  invested  in  improving  the  model  for  a  given  enclosure  and  exposure  scenario. 
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page  1 

PROTECTION  FACTORS  FOR  SOLID  AEROSOL  PENETRATION  INTO  SEALED  ENCLOSURES 

PFACDOC.MCD  Version  19  Feb 

Current  test  ID;  1221Tlx 


PARAMETER  VALUES: 


particle  density 
indices 


p  *  2 .7  g/cni3 


ORIGIN  s 1 


a  B  1  .  .  11 
X  s  12  .  .  20 


j  s  1  ..  5 
O  a  1  ,3  .  .  20 
i  t  1  .  .  20 


sizes  di££, 
microns !  press 

mjj  dpWGj 


Exp .  Exp . 
sizes  decay 

Nj  exBj 


sizes 

dj  a  m^*  10'* 

cm 

.1 

2 

.2“ 

3 

gravity  acceleration 

g  B  980 

cm/sec2 

.3 

±_ 

.4 

5 

air  viscosity 

h  a  1.83'10'^ 

g/cm/s 

.5 

6 

.6 

7 

Boltzmann  constant 

k  a  1.38-10"^* 

erg/K 

.7 

0 

.8 

9 

Temperature 

Temp  a  288 

deg  K 

.9 

10 

1 

flow  exponent 

n  a  .5 

1.5 

flow  coefficient  Cd ■  1 


Turbulent  dissipation  parameters: ke B 1000  [l/sec]  ne  b 2 
empirical  multiplier:  mult  =  ,189 

Dimensions  of  rectangular  enclosure: 

Height:  Length:  Width; 

LI  i  117  L2  =  213  L3  =  167  [cm) 


aperture  depth:  height  width:  perimeter:  area: 

1  B  4  h  B  0 .015  w  ■  80  P  B  2*  (w  ■f  h)  cm  Ar  ■  wh  cm2 

smallest  dimension:  bsh 


Cloud  duration: 


Ts  .08333 


Exposure  time  tsO.B  [hours] 


Lognormal  spectrum  parameters: 


Sigma:  Geom.  mean: 

sigGB  2.284  meanGBi,2 
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PRELIMINARY  CALCULATIONS: 


-2- 


air  flowrate: 


Or.j  -J  406-Cd'Ar'  cm3/s 


Cunningham'S  correction: 

Pd,  s  1  +  •exp(-7.44'10'*‘di) 

d,  ' 


d, 


Or. 

Vr,  s  — i  cm/sec 
^  Ar 


Terminal  velocities: 

Diffusion  coefficients: 
Settling  parai.ieter: 


Db,.J^®n!P-.Pdi 

‘  3-K-Tl-di 


XI. 


R'Ut, 


X 

n* 


Reynold '  s  n'-.Tibera : 


Lognormal  distribution: 


lognorn,  ■ 


Rn^  ab-Vrj-1.2 


10 


.  3 


1 

. — I..,.,.  I, 

■V(2 •  Jt)  •  In  ( sigG)  -mj 


•  exp 


In 


"'i  \\\* 


meanO 


In(BlgO)  / 
-  2 


Rn, 


40 


PENETRATION  TRANSPORT  -  DEPOSITION  BY  MOLECULAR  DIFFUSION: 


Under  laminar  flow  conditions  the  deposition  of  particles  is  controlled 
by  the  diffusion  coefficients,  the  length  of  the  channel  and  the  flow 
through  the  channel.  The  controlling  parameter  is  then  defined  as: 


P 


i  .i 


■  X'Db 


For  transport  through  a  rectanqular  channel: 
foSaj^j  *  .819-exp(-3.65-p^^j)  +  .097-exp(- 22  .  3*P^ 


faSbj  B  .032-exp (- 57- +  .027-exp(- 123-Pj_j)  +  . 025- exp(-  750- P^  _ 
foBDR^  _  j  e  faSa^  +  fa5b^  ^^ 

If  P<0.1  use: 

fa6DRj_jEl  -  2.56-(Pj,  +  l-2-Pv,.l  +  •  ’  (Pi .  j) 
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3- 


cm-l 


PENETRATIOH  TRANSPORT  >  TURBULENT  IMPACTION i 


Deposition  will  be  determined  by  turbulent  eddy  transport  and  the  inertial 
properties  o£  tho  particles.  The  deposition  velocity  must  be  derived 
from  considerations  of  Reynold's  number  and  surface  roughness.  But  using 
the  terminal  settling  velocity,  Ut,  for  now  leads  to  minimum  wall  loss. 


microns 
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-4- 


PENETRATION  TRANSPORT  -  GRAVITATIONAL  SETTLING: 


Settling  is  dependent  on  leak  rate,  particle  density  and  diameter,  and 
channel  dimensi.ons , 

Using  Thomas'  egn.  for  gravity  settling  in  a  rectangular  channel: 


fttGR,  ,  s  1  -  l*Uti*  — 
i-J  ^  Qr^ 


‘oTT 

0.3 

o.s 

0.7 


0^9 

1.5 


9 


microns 


Selected  transport  fractions; 


£aaR|^  fotGR^o.^) 


1 

0.993 

1 

0.995 

1 

0.995 

1 

0.996 

1 

0.996 

fqQRny 

0.562 

0.642 

0.69 

0.723 

0.747 


TOTAL  TRANSPORT  FRACTION: 


The  total  transport  fraction  is  defined  as  the  ratio  of  the  total 
number  of  particles  exiting  the  aperture  to  the  number  entering. 
The  total  number  exiting  will  be,  in  terms  of  the  above  transport 
fractions,  the  number  entering  less  the  loss  due  to  each  process; 


Nout  ■=  Nin  -  (1  -  faD)  -Nin  -  (1  -  foT)  -N^  -  (1  -  fOG)  -Nt^a 

The  total  transport  fraction,  Nout/Nin,  reduces  to: 
f  :=  faD  +  faT  +•  faG  -  2o 

For  P>0.1  use  fR^  j  a  faSDRj  +  faTR^  ^  +  faOR^^^  -  2o 

For  P<0.1  urie  fR^  »  faLDRj  4-  foTR^  ^^  f  faGR^  ^^  -  2 


Appendix  B 


65 


-5- 


TOTAL  PENETRATION  TRANSPORT  MATRIX; 


0.999 

1 

1 

1 

1 

0,999 

0.999 

0.999 

0.999 

0.999 

0.997 

0.998 

0.998 

0.998 

0.999 

0.996 

0.997 

0.997 

0.997 

0.998 

0,994 

0.995 

0.996 

0.996 

0.997 

0,992 

0.994 

0.994 

0.995 

0.995 

0.99 

0.992 

0.993 

0.993 

0.994 

0.9  87 

0.989 

0.991 

0.992 

0.992 

0.984 

0.987 

0.989 

0.99 

0.991 

0,981 

0.984 

0.986 

0.988 

0.989 

0.959 

0.967 

0.971 

0.974 

0.976 

0.93 

0.943 

0.951 

0.956 

0.96 

0.852 

0.  878 

0.894 

0.905 

0.913 

0.747 

0.791 

0.818 

0.837 

0.85 

0.621 

0. 685 

0.725 

0.752 

0.773 

0.477 

0.S63 

0.617 

0.654 

0.682 

0.319 

0.428 

0.497 

0.544 

0.58 

0.152 

0.283 

0.366 

0.425 

0.469 

■0.022 

0.131 

0.228 

0.297 

0.35 

■0,2 

-0.027 

0.084 

0.164 

0.224 

dpWGj 

1 

1.5 
2 

2.5 
3 


microns 


(Negative  values  reflect  exaggerated  loss  of  largest  particles) 


microns 


(Differential  pressure  increases  with  higher  curves) 


Appendix  B 


66 


STIRRED  SETTLING  OP  PARTICLES: 


-6 


N.B.-  Turbulent  eddy  diffusion  is  defined  by: 
Particle  decay  rates  for  rectangular  chamber  [1 


De  =  ke- Y"*a 
/sec]  . 

ke  =  1*10^ 
ne  >  2.038 


microns 

Calculate  theoretical  lognormal-weighted  decay  rate: 

^  lognontii 

^lognormi-m7“  InB  =  6.922-10“*  1/sec 

i 


Calculate  experimental  lognormal -weighted  decay  rate: 
(Experimental  particle  decay  rates  derived  from  1204T7  test  dita) 


wexB  * 


lognorm^'exBi’mj 

i _ 

y lognorm^‘m^ 
i 


wcxB  =  6.692*10 


1/sec 
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-7 


PROTECTION  FACTORS  FOR  THE  EXPOSURE  SCENARIO: 


Air  infiltration  rate: 


Qr.-3600 

Rj  ■  7— ’tt-  /:■  [1/hour] 

LI  *  L2  *  L3 


Total  particulate  removal:  +  Bi*3600  [l/hour] 

Display  air  flowrates  and  selected  removal  rates: 


mmWQ  [m3 /hr]  [1/min] 


[l/hr]  tl/hr]  [1/hr]  [l/hr] 


0.1  1.0 


10  .microns 


Calculate  Protection  Factors: 
Vapor: 


Pv 


i 


Rj-T 

Rj‘T  +  exp(;^- R^‘t)  •  (^1  -  exp(^Rj-T)]) 


Particulate: 


Pp 


1  .J 


[RPi,^-T  +  exp(-Rpi,j-t).(l  -  exp(RPi,j-T))  ] 


RP 


i  .J 


Eliminate  negative  values  resulting  from  exaggerated  filtration: 


PPi.j  “  0  ,Ppj  _  ,PPi,j') 
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-8- 


THEORETTCAL  PROTECTION  FACTORS:  l-i2lTl  talc/nozzle 
dpWG^ 


1. 


1 . 5 


2.5 


p  =  2.7 


1  = 
h  = 
w  = 


4 

0.015 

80 


[mmWG) 


Particles i 


Vitpor  I 


Pv 


/  5'n 

4 . 7 
4.2 

3.8 

3.5/ 


[cm] 

Icifl] 

[cm] 


Pp 


Cd  =  1 

T 

:=  0.083  [hr] 

n  =  0. 

5  t. 

=  0.5 

[hr] 

ke  =  1 

•10^  ne  =  2 . 

038 

17.4 

14.3 

12.5 

11.3 

10.4 

11.6 

9.6 

8.4 

7.6 

7 

9.8 

8.2 

7.2 

6.5 

6 

9 

7.5 

6.6 

5.9 

5.5 

8.5 

7.1 

6.2 

5.6 

5.2 

8.2 

6.8 

6 

5.4 

5 

8 

6.S 

5.8 

5.3 

4.9 

7.8 

6.5 

5.7 

5.1 

4.8 

7.7 

6.4 

5.6 

5.1 

4.7 

7.6 

6.3 

5.5 

5 

4.6 

7.7 

6.4 

5.6 

5.1 

4.7 

8.5 

7 

6 . 1 

5.5 

5 

11.6 

9.3 

8 

7.2 

6.6 

17.5 

13.7 

11.6 

10.2 

9.3 

28.5 

21.3 

17.6 

15.3 

13.7 

49. 9 

34.8 

27.7 

23.5 

20.7 

98.2 

60.1 

45.2 

37.1 

31.9 

264.9 

116.4 

78.4 

60.7 

50.3 

264.9 

314.9 

157.1 

108.2 

84.1 

264.9 

314.9 

519.3 

239.8 

160 

microns 
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WEIGHTED  PROTECTION  FACTORS  FOR  THE  GIVEN  SIZE  SPECTRUM: 


pressure  [nunWG] 
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-10- 


WEIGHTED  PROTECTION  FACTORS  FOR  LOGNORMAL  DISTRIBUTION: 


Distribution 

parameters 

definition  (dp) 

P 

sigQ  s  2.284 

meanG  s  i .  2 

Mode 

-1 

Geometric  mean  or  mode 

0 

chosen  p; 

Arithmetic  a«an 

0.5 

Diameter  of  average  area 

1 

p  B  -  1 

Diameter  of  average  mass 

1.5 

Surface  median  diameter 

2 

dp  s  meanG-exp 

C.P‘  In  ( sigG)  V 

Surface  mean  diameter 

2.5 

Volume  (mass)  median  diameter 

3 

Volume  (mass)  mean  diameter 

3.5 

dp  =•  0.607 

microns 

nig  lognonUg 


o.i  10 


0 

.052 

0 

.394 

0 

.551 

^0 

.558 

.505 

■j 

.31 

0 

.087 

0 

.022 

0 

.007 

n 

.003 

i H  1  .  .  16  ftruncace  at  m  <  7  microns) 


^  lognomii-Ppii  -nii 

PPlllj  *  ■  —  . 

Vloanormj-m^ 

i 


Vapor;  Particles; 


pririj 


dpwa^ 
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PROTECTIOfl  FACTORS  USING  EMPIRICAL  SETTLING  DATA; 

Settling  rate  corrected  relative  to  original  datat 
exRpj^j^Rj  -t-  exB^'3600*mult 


0.1  ntj  10 


microns 


exPPi , j  e  — 


exRpj ,  j*T 


exRp^^j'T  +  exp(]- exRp^  1)‘ (^1  -  exp(jaxRp^  'ij. 

exRp^ , j 


Eliminate  negative  values  from  results; 

exPpi  ^  »  i £  (Re  (exPPi  0  ,  exPPi  _  i ,  j  ,  exPp^  _ 


Calculate  experimental  lognonnal  weighted  protection  factors: 


exPFlnj 


^  lognomii-exPPu^j, 


i 
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APPENDIX  C 

Utility  Data  Analysis  Programs 


C.l, 


Data  acquisition  from  the  LAS-X  Aerosol  Spectrometer  is  performed  by  a 
dedicated  ltP-8SA  personal  computer  and  analyst  program  (Soderholm,  1984).  Since  data 
analysis  will  eventually  be  performed  on  an  IBM-compatible  machine,  it  was  decided  to  transfer 
the  data  to  a  Mathcad  program  which  would  provide  the  functions  of  background  subtraction, 
dilution  correction  and  spectrum  display.  Later  versions  could  be  adapted  for  automated  data 
acquisition  and  alternative  spectral  characteristics,  such  as  lognormal  fitting  and  mass-related 
fuctions.  The  LAS-X  provides  four  size  ranges,  each  comprising  1 5  individual  bins,  and  an 
oversize  bin  for  each  range.  The  size  ranges  are  as  follows. 


Range  3:  0.120  -  0.210  micron 
Range  2:  0.170-  0.590  " 

Rangel:  0.300-1.800  " 

Range  0:  1.500-  7.500  " 


interval  =  0.006  micron 
''  0.03 

•'  0.10  " 

"  0.40 


The  size  intervals  of  the  ranges  differ  widely  one  from  tiie  other,  hence  it  is  difficult 
to  visually  grasp  spectra  wliich  span  several  ranges.  A  utility  program  was  writf  'n  in  Mathcad 
language  to  facilitate  LAS-X  data  analysis  by  compressing  the  background-  and  dilution-corrected 
counts  from  each  range  into  new  intervals  which  are  easier  to  visualize.  Results  are  displayed  in  a 
new  .c.trum  comprising  42  bins  instead  of  60.  The  compression  ratios  of  the  smaller  size  bins 
can  be  easily  adjusted  to  account  for  any  spectra  of  interest  in  the  future.  The  following  program, 
named  LASXDOC.MCD,  is  provided  with  minimum  comments  in  the  body  of  the  text. 

Additional  comments  are  provided  here 


Pakc  I: 


The  program  is  initialized  with  identifying  numbers  for  the  test.  TID  and  BID  are  the 
sampling  times  automatically  printed  on  the  LAS-X  outputs  with  each  run.  The  dilution  factor 
x20  is  included  if  the  TSI  Diluter  was  used  in  the  run.  If  not,  the  number  1  is  written.  The  indices 
i  through  q  determine  the  subdivisions  -^f  the  ranges  that  are  to  be  compressed.  They  can  be 
changed  to  tailor  the  program  for  specific  solid  aerosols 

INPUT  BACKGROf  IND  DATA  Of  THE  FOUR  SIZE  PANGES 

The  dxj  vectors  are  the  fixed  sizes  of  each  bin  in  the  four  ranges.  The  bkgdxj  vectors 
are  the  raw  data  counts  H  orn  each  bin  printed  out  by  the  HP-85A  program.  They  are  replaced  by 
the  new  date  at  each  run  If  the  same  background  run  is  used  repeatedly  for  a  scries  of  runs,  they 
do  not  have  to  be  altered 

C'ORKECT  FOR  x2.0  DILUTION  AND  EFFICIENCY  IF  APPLICABLE 


74 


This  segment  corrects  each  background  range  for  the  dilution  factor  if  dil=20.  The 
highest  size  range,  bkgdO},  includes  a  correction  factor  for  the  efficiency  of  the  diluter  at  these 
sizes.  The  manufacturer's  efficiency  curve  was  previously  fit  to  a  polynomial  function  using  the 
ASYSTANT  curve  fitting  option, 

Page  2: 

CALCULATE  COMPRESSED  BINS  FOR  RANGE  3; 

Range  3  is  compressed  into  three  bins,  each  including  five  original  bins  as  determined 
by  the  indices  j,  k  and  1.  The  interval  size  is  now  0.018  micron. 

DEFINE  NEW  BIN  SIZES  FOR  THE  COMPRESSED  RANGE  3 : 

These  constants  are  the  new  mid-point  sizes  redefined  for  the  compressed  range  3. 

CALCULATE  COMPRESSED  BINS  FOR  THE  SECOND  HALF  OF  RANGE  1 : 

This  segment  compresses  the  range  1.0-1. 7  microns  to  provide  a  smoother  transition 
to  the  highest  range  0. 

DEFINE  NEW  SIZE  BINS  FOR  THE  COMPRESSED  RANGE  1 : 

These  constants  are  the  new  mid-points  redefined  for  the  compressed  range  1 . 

PLOT  THE  COMPRESSED  BACKGROUND  DATA  COUNTS: 

The  plot  shows  the  smoothed  background  spectrum  after  compression  of  range  3. 
Ambient  background  at  the  test  site  is  always  below  0.5  micron  so  that  range  0  is  not  affected. 

Page  3 

INPUT  EXPERIMENTAL  DATA: 

The  four  dataxj  vectors  are  the  raw  count  data  of  the  test  run  itself  They  correspond 
to  the  size  ranges  appearing  on  page  1 ,  They  are  replaced  by  new  data  vectors  for  each  run 
whether  or  not  the  background  data  had  to  be  changed.  Recall  that  the  test  identification  number 
and  the  aerosol  name  are  to  be  updated  for  each  run.  It  is  useful  to  also  change  the  file  name  each 
time  to  indicate  the  nature  of  the  run  and  to  save  results  individually. 

CORRECT  DATA  FOR  DILUTION  FACTOR  AND  EFFICIENCY  IF  APPLICABLE: 

As  for  the  background  data,  the  test  data  is  adjusted  for  dilution,  if  needed. 


Appendix  C 


75 


CALCULATE  COMPRESSED  BINS  FOR  RANGE  3/  SECOND  HALF  OF  RANGE  1 : 

As  for  the  background  data,  the  test  data  are  compressed  in  a  like  manner.  The 
sizes  of  the  new  bins  were  calculated  previously. 

Page  4: 

PLOT  COMPRESSED  DATA  COUNTS: 

The  test  data  are  displayed  with  the  new  size  bins.  Note  that  there  is  still  overlap 
between  ranges,  but  the  plot  is  smoother  than  the  original  uncompressed  version. 

CALCUI.ATE  NET  COUNTS  ABOVE  BACKGROUND/PLOT  NET  COUNTS; 

Background  counts  in  the  compressed  ranges  are  subtracted  from  the  data  counts  to 
provide  net  counts  in  each  range  and  plotted. 

Page  S: 

The  compressed  net  counts  in  the  new  bins  are  tabulated  alongside  the  newly 

defined  sizes. 
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Page  1 

LASXDOC.MCD 

CALCULATION  OF  NET  PARTICLE  COUNTS  FROM  LAS-X  SPECTROMETER  DATA 

version  23  Sept  92 

This  program  accepts  raw  size  count  data  from  the  PMS  LAS-X  Laser  Aerosol  Spectrometer 
and  com.'erts  them  to  a  compressed  background-  and  dilution-corrected  spectrum. 

Aerosol:  Ti02Tiona  Test  identification  TID  =  150540 

Input  dilution  factordll  20  Background  identification  BID  =  131130 


i:=  l„15 
n  ;=  1..6 


j:=l„5 
0  :=  7, .9 


k  ;=  6..  10 

p  :=  10..  12 


1-  11..  15 
q  :=  13..  15 


INPUT  BACKGROUND  DATA  OF  THE  FOUR  SIZE  RANGES: 


d3; 


bkgd3j 


d2i  = 


,120 

,126 

,132 

,138 

,144 

.150 

.156 

.162 

.168 

.174 

,180 

.186 

.192 

.198 

.204 

125 


120 


112 


129 


109 


119 


136 


121 


116 


84 


97 


89 


70 


65 


66 


17 


20 


23 


26 


29 


32 

35' 


38 


41 


44 


47 


50 


53 


56 


59 


bkgd2.  = 


dL 


258 


168 


64 


36 


18 

11 


.5 


,6 


.8 


.9 


1.0 


1.1 


1.2 


1.3 


1.4 


1.5 


1^ 

1.7 


m:=  1..3 
ORIGIN  il 


bkgdlj  :=  do.  bkgdO.  := 


12.3" 


P.7 


3.1 

3,5 


3,9 


1.5 


K.3 


W.7 


5.1 


5.5 


5.9 


16.3 


6.7 


7.1 


0 

0 

0 

0 

0 

0 

I 

0 

p 

0 

0 

p 

0 


micron  counts  micron  counts  micron  counts  micron  counts 

CORRECT  FOR  x20  DILUTION  AND  EFFICIENCY  EF  APPLICABLE: 

bkgd3j  :=  if^dil  -  20,bkgd3i,bkgd3j  dil^  bkgd2j  =  if^dil  -  20,bkgd2j,bkgd2-dil^ 


bkgdlj  :=  if^dil  -  20, bkgd  1  j, bkgdlj' dil^ 

bkgdOj  ifjdil  -  20,bkgd0j,bkgd0j  dil'[  1,06  -  .024- dO^  -  .002' (dOj)^]] 
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CALCUIATE  COMPRESSED  BINS  FOR  RANGE  3: 

bkgdSxj  ;=  ^bkgd3j  bkgdSxj  =  ^bkgd3^  bkgd3x3  =  |^bkgd3| 

j  k  1 

DEFINE  NEW  BIN  SIZES  FOR  THE  COMPRESSED  RANGE  3: 
d3x,;=.132  dSxj  :=  .162  dSxj  :=  .192 

CALCULATT  COMPRESSED  BINS  FOR  THE  SECOND  HALF  OF  RANGE  1: 

bkgdlx^  :=y^bkgdlq  bkgdlXj  =  ^bkgdlp  bkgdlxj  =^bkgdlj 

op  q 

DEFINE  NEW  BIN  SIZES  FOR  THE  COMPRESSED  RANGE  1 : 
dlXj  =  1,0  dlX2  :=  1.3  dlXj  :=  1.6 

PLOT  THE  COMPRESSED  BACKGROUND  DATA  COUNTS: 


13000 

bkgdSXjjj 

JL 

bkgd2. 

bkgdl^ 

1. 

bkgdlXj^ 

JL 

bkgdO. 

ji. 

0 


BID -131130 


rs 

■ 

If 

■ 

1 

1 

1 

II 

■d3x^.d2j.dl^.dlx^:d0i  10 


microns 
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INPLT  EXPERIMENTAL  DATA: 

data3|:=  data2j:=  datalj-  dataOj  = 


557 

369 

285 

308 

269 

270 
282 
294 
218 
292 
205 
226 
236 
196 
205 


1855 

1052 

769 

755 

654 

674 

568 

558 

471 

518 

507 

611 

935 

kiO 

2076 

3375 

3567 

9237 

10969 

6386 

4441 

3261 

2632 

1982 

1556 

1306 

959 

691 

558 

441 

3379 

1454 

638 

332 

176 

109 

43 

31 

18 

10 

12 

2 

2 

3 

2 

counts  counts  counts 


counts 
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CORRECT  DATA  FOR  DILUTION  FACTOR  AND  EFFICIENCY  IF  APPLICABLE: 
data3|  =  if^dil  -  20,data3.,data3|'dil^  data2|  if^dil  -  20,data2.,data2.'dil^ 


dntalj  ;=  if^dil  -  20,datnl.,dalalj  dil^ 

dataOj  :=  if|dil  -  20,data0j,data0j  dil'  1.06  -  .024  d0j  -  ,002-^d0j^^ 


CALCULATE  COMPRESSED  BINS  FOR  RANGE  3: 

data3xj  ;=  ^data3.  dataSXj  :=  ^data3jj  dataSxj  :=  ^data3j 

j  k  1 

CALCULATE  COMPRESSED  BINS  FOR  SECOND  HALF  OF  RANGE  1 : 


datalxj  :=  ^dataljj  datalxj  -  ^datalp 

0  p 


datalXj 


^datal, 
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PLOT  COMPRESSED  DATA  COUNTS: 


Page  4 


1  d3x^.d2i.«il,.dlx^.d0,  10 


microns 

CALCULATE  NET  COUNTS  ABOVE  BACKGROUND: 
netSo,  datajx^j  -  bkgdSx^j^  net2j  :=  data2j  -  bkgd2| 

netljj  :=  dataljj  -  bkgdl^  "  bkgdO. 

PLOT  NET  COUNTS: 

TID  =  150540  Ti02Tiona 


microns 
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Page  5 


DISPLAY  COMPRESSED  NET  COUNTS; 


d3x_ 

m 

net3^ 

0.132 

23860 

0.162 

15600 

0.192 

13620 

micron 

counts 

d2.  net2. 


micron  counts 


din 

netl„ 

0.3 

67280 

0.4 

71280 

0.5 

184720 

0.6 

219380 

0.7 

127720 

0.8 

88820 

micron 

counts 

dlx„  netlx_ 
m  m 


1 

1.3 

157480 

76400 

1.6 

33800 

microns  counts 


dOj 

netOj 

1.5 

68898 

1.9 

29289 

2.3 

12686 

2.7 

6511 

3.1 

3402 

3,5 

2074 

3.9 

805 

4.3 

570 

4.7 

325 

5,1 

177 

5.5 

”208 

5.9 

34 

6.3 

33 

6.7 

49 

7.1 

32 

microns 

counts 
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C.2.  Exponential  Decay  Rate  Calculation 

The  estimation  of  exponential  decay  rates,  whether  A'om  tracer  gas  concentration  or 
solid  aerosol  settling  in  a  stirred  enclosure,  is  done  by  least  squares  fit  of  the  logarithm  of  the 
concentration  data  against  time.  A  utility  program  written  in  Mathcad,  LASXRATE.MCD  shown 
on  the  next  page,  serves  all  applications  of  this  nature. 

The  example  shown  is  for  the  settling  of  particles  from  a  certain  size  bin  of  LAS-X 
data.  First,  the  number  of  data  points,  N,  is  specified.  A  time  vector,  time|c,  is  inputted 
corresponding  to  the  measurements,  count);.  Background  counts,  if  constant  over  the 
measurement,  are  inputted  as  bkgd);.  If  the  measurement  is  of  long  duration,  such  that  the 
background  is  changing  as  well,  a  Unction  can  be  inputted  instead. 

Net  counts  are  calculated  and  normalized  by  the  first  value  for  simplicity.  Mathcad 
functions  INTERCEPT  and  SLOPE  then  calculate  values  for  the  zero  intercept,  exp(u),  and  the 
slope,  V.  The  correlation  coefficient  is  calculated  by  Mathcad  function  CORR.  The  decay  rate  in 
units  of  sec~^>  in  this  case.  The  normalized  data  points  and  the  fitted  function  are  plotted 
together. 

The  utility  program  may  be  used  for  any  exponentially  decaying  data  set  as  long  as 
consistent  time  units  are  maintained. 
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Data  index: 

timejj:= 

T 

4' 

L 

\2_ 

16_ 

W 

28 

32 

36 

40 

11 

48 

52 


EXPONENTIAL  DECAY  RATE  CA1,CULATI0N 
LASXRATE.MCD 
data;  1204t7  bin  3.90-4.30  micron 

N  :=  14  points  lc:=l..N  Background;  bkgd|j  =  0  counts 
countjj  :=  Calculate  net  counts  and  normalize: 

countjj  -  bkgdjj 

cv  := - InC,  :=  In  &. 

countj  ^ 

Calculate  best  lit  to  exponential  function: 

u  =  intercept  (time.lnC)  u=0,12 

exp(u)- 1.128 

V  :=  slope(time,lnC)  v  =-0.0573 

Calculate  correlation  coefficient; 

corr(time,InC)  =-0.95857 
Calculate  decay  rate  in  1/seconds: 

rate  =  —  rate  =  9.55*  10  ^  sec-l 

60 

1 


'Tc 

0 

cxp^uf 


0 

minutes 


3l] 

ir 


US 
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C.3.  Theoretical  Settling  Rate  Calculation 

In  order  to  study  the  theoretical  model  of  particle  settling  in  a  stirred  enclosure 
independently,  the  utility  programs  SETTLHivEDOC.MCD  and  SETTLEN.DOC  were  compiled, 
These  are  essentially  the  same  as  the  submodel  that  is  integrated  into  the  protection  factor  model 
program  but  used  separately  to  allow  parameter  sensitivity  tests  with  the  turbulent  dissipation 
intensity,  kg,  and  the  exponent,  n.  an  an  independent  parameters.  Any  other  parameter  can  be 
tested  if  the  appropriate  changes  in  the  indices  are  made.  Sample  runs  of  the  two  programs 
folow. 


The  dimensions  of  the  rectangular  enclosure  are  those  of  the  test  exposure  chamber. 

The  remaining  parameter  values  given  are  typical.  The  particle  settling  rates  are  calculated  in  a 
two-dimensional  matrix  as  a  function  of  size  and  the  chosen  parameter  vector. 

Results  are  plotted  on  Page  2  of  each  run,  Terminal  velocities  and  diflusion 
coefficients  are  plotted  for  reference,  In  the  first  run,  the  settling  rate  matrix  for  varying  turbulent 
intensity,  k^,  (with  constant  exponent  n-2.0)  is  plotted  below  with  the  five  runs  clearly 
delineated.  Note  that  below  1  micron  the  settling  rate  varies  by  a  factor  of  about  3  for  a  change  in 
kg  of  10.  Above  3  microns  the  settling  rate  is  insensitive  to  the  turbulent  dissipation  rate.  The 
second  run  shows  the  effect  of  changing  n  with  a  fixed  value  of  kg,  A  strong  sensitivity  to  the 
exponent  n  is  evident  as  the  settling  rate  below  1  micron  varies  by  about  an  order  of  magnitude 
for  a  change  in  n  from  2  to  2.8,  the  range  in  which  published  values  have  been  noted.  Again,  the 
settling  rate  above  3  microns  is  insensitive  to  this  parameter.  In  this  region,  gravitational  settling 
predominates.  The  particle  size  with  the  minimum  settling  rate  is  very  sensitive  to  the  parameter. 
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MODEL  EQUATIONS  FOR  STIRRED  SETTLING  OF  PARTICLES 
(SETTLEKE.MCD) 


Dimensions  of  rectangular  chamber: 

indices: 

height- 

length- 

width- 

isl..32  pb24 

L1*117 

L2s213 

L3 ■  167  cm 

jel..p  Np..32 

Parameter  values: 

P*1 

g/cm3 

Ti298 

degK 

gc980 

cm/sec2 

n.2.0 

TIB  1, 83- 10“* 
kil.3810''® 

g/cm/s 

erg/K 

dj«mj'  10'* 

(conversion  to  cm,) 

Cunningham's  corrections: 


,  1,76-10'®  5, 6. 10'^ 

Fd|S  1 4- - - - +■ 


d: 


di 


exp 


(-7,44-10'^'di 


Terminal  velocities: 


ut|.-eA.(di)4F(j| 
'  18.11  ^  */  ‘ 


Diffusion  coefQcients: 


k-T 

Uh,^-JL± — Fd. 
3  X  Ti  dj 


Settling  Parameter: 


71- Ut; 


"■“"S 


kc,.(Dbi) 


n-  1 


Particle  settling  rates  for  rectangular  chamber  [LsecJ: 


XI.  q  Il2  L3/ 


Ut: 


Lltanhf — 

\  \  2  II 


qsl..5 


100 

250 

500 

750_ 

Tooo 


sec-1 


sizes: 


.05 

.06 

.07 

.08 

.09 

.1 

.125 

.15 

.175 

.2 

.25 

.3 

.4 

‘J" 
.6 
.7  ' 
.8 
.9  ' 
1 

1.2 

1.4 
1.6 
1.8 
2 

2.5 

3 

3.5' 

4 

5 

6 
8 
10 


microns 
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0.01  0,1  1  10 


microns 


Settling  Rntes  fl/sec] 


0.01  0.1  1  10 


100 

250 

500 

750 

MO^ 


m- 


microns 
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MODEL  EQUATIONS  FOR  STIRRED  SETTLING  OF  PARTICLES 
(SETTLEN.MCD) 


Dimensions  of  rectangular  chamber: 
height-  length-  width- 

Ll«117  L2«213  L3il67  cm 


indices; 

i»1..32  p»24  q-1..5 

l*p.,32 


Parameter  values; 

pil 

g/cm3 

T-298  degK 

gB980 

cm/sec2 

ke>1000 

ti«1.83'10*'‘ , 

g/cm/s 

dj  ■  m.- 10'^  (conversion  to  cm.) 

lcsl.38'10‘^ 

erg/K 

Cunningham's  corrections; 


2 

U 

U 

2.8 


_  ,  1.76-10‘*  5,6' 10'®  .^.4. 

Fd|«l+ - : - + - : - exp(-7.44-10  -d, 


d, 


d, 


Terminal  velocities: 


Uti 


P'8  _f.\2 


la-ri 


(-“i)  -Mi 


Diffusion  coefficients; 


k'T 

Db  .-Li — Fdj 
Sjc-Ti'dj 


Settling  Parameter; 


Tl-Uti 


Particle  settling  rates  for  rectangular  chamber  [1/sec]: 


Pi 


2Ut, 


Xl^  ,  \U.  u 


Ul. 


Ll-tanh — 
2 


Page  1 


sizes: 


.05 

.06 

.07 

.08 

.09 

.1 

.125 

.15 

.175 

.2 

.25 

■3 

.4 

.5 

.6 

.7 

.8 

.9 

1 

1.2 

1.4 
1.6 
1.8 
2 

2.5 

3 

3.5 

4 

5 

6 
8 
10 


microns 
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Page  2 

Terminal  Velocities  and  Diffusion  Coefficients 


microns 


Settling  Rates  [1/sec] 


2 

11 

2.4 

2.6 

2.8 


tn. 


microns 
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Sample  results  of  the  attempts  to  optimize  the  two  parameters  of  the  turbulent 
dissipation  coefficient  were  given  in  the  text.  The  optimization  was  performed  using  the  Asystant 
scientific  software.  Aside  ftom  the  option  to  define  user  functions,  the  program  is  entirely  menu- 
driven.  In  order  to  insert  the  settling  rate  model  for  curve-fitting  against  the  data,  the  equations 
were  defined  as  a  five-stage  nested  user  fliction.  On  the  following  two  pages  the  screen  prints 
showing  the  user  function  definitions,  parameters  and  variables  are  shown.  Fourteen  experimental 
data  points  were  inserted  in  variable  Y  against  their  sizes,  inserted  in  variable  Z.  Prior  to  the 
optimization,  Z  and  Y  are  placed  on  the  stack  followed  by  the  final  function,  F7.  The  result  of  the 
calculation  was  stored  in  variable  B,  labelled  beta. 

The  nested  function  definitions  are  as  follows: 

FI  -  Cunningham's  correction 

F2  -  Terminal  velocity 

F3  -  Diffusion  coefficient 

F4  -  Settling  rate  parameter  (XI) 

F7  -  Settling  rates 

At  the  bottom  of  the  second  page  is  the  screen  for  chosing  the  fit  of  the  experimental 
data  (variable  Y  with  the  sizes  Z),  against  the  curve  calculated  by  function  F7  with  two 
independent  parameters,  A  and  B,  whose  initial  values  are  given  in  the  leff-hand  column. 

The  third  page  shows  screens  for  three  results  given  in  the  text.  Note  that  for  widely 
varying  kg  parameter  values  (1000<A<  10000),  similar  values  of  n  (2.038<B<2.301)  yield 
essentially  identical,  and  rather  unsatisfactory,  fits  to  the  experimental  data.  This  is  despite  the 
fact  that  the  least  squares  fit  correlation  coefficient,  R^,  is  excellent  in  all  cases 
(0.9335<r2<0.9361).  The  model  as  posed  is  ill-conditioned  for  fit  through  these  two  parameters 
against  the  measured  data. 

Although  the  model,  in  this  case,  does  not  result  in  satisfactory  results,  the  use  of 
Asystant  as  a  tool  in  quantitative  sensitivity  analysis  is  demonstrated  as  a  useful  time-saver  for 
whatever  models  may  be  developed  in  the  future. 
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— — *  Function  tSefinitions  «  — — 

Press  Esc  when  done 
FI  1>  \  l+1.76E-5/Z+5.6E-6**xp{-7.44E-4*Z)/Z 
2> 

3> 

4> 

5> 

Use  PgUp  and  PgDn  to  edit  other  functions 
I  —■  Calculator  Functions 


Parameters 


store 

drop 

Bsin 

acos 

Btan 

in 


pedlt 

swap 

alnh 

cosh 

tanh 

lOexp 


Main  Menu  Options 


Help 

File  Proc 
Wave  Proc 
Wave  Gen 


Graphics 
File  I/O 
Curve  Fit 
Stats 


status  print 

roll  pi 

asinh  Inv 

acosh  neg 

atanh  abs 

log  sgrt 


Functions 
Polys 
Diff  Eqs 
Save/Exlt 


A-  1000 
B-  2.S28781E0 
C-  0 

D-  O.OOOOOOE-l 
E-  O.OOOOOOE'l 
F-  O.OOOOOOE-l 
G-  O.OOOOOOE-l 
H-  O.OOOOOOE-l 
I-  O.OOOOOOE-l 

■  Variables 


rho 

R- 

2.700000EO 

beta 

S- 

14  DP-REAL 

eta 

1.830000E-4 

T 

u- 

298 

LI 

V- 

117 

L2 

w- 

213 

L3 

X- 

167 

data 

y- 

14  DP-REAL 

sizes 

z- 

14  REAL 

•mmmmm  Functlon  Definitions 

Press  Esc  when  done 
1>  \  S*9aO*2‘2*Fl/ (18*T) 


Use  PgUo  and  PgDn  to  edit  other  functions 
Calculator  Functions  . . 


next 

store 

psdit 

status 

dup 

sin 

drop 

Bsln 

swap 

Sinn 

roll 

asinh 

cos 

acos 

cosh 

acosh 

tan 

atan 

tanh 

atanh 

exp 

In 

lOexp 

log 

Main  Menu  Options 


Kelp 

File  Proc 
Wave  Proc 
Wave  Gen 


Graphics 
Fils  I/O 
Curve  Fit 
Stats 


aqrt 


Functions 
Polya 
Diff  Eqe 
Save/Exit 


rho 

beta 

eta 

T 

LI 

L2 

L3 

data 

sizes 


■  Parameters  mmtm 

A-  1000 
B-  2.528781E0 
C-  0 

D-  O.OOOOOOE-l 
E-  O.OOOOOOE-l 
F-  O.OOOOOOE-l 
G-  O.OOOOOOE-l 
H«  O.OOOOOOE-l 
!■  O.OOOOOOE-l 

■  Variables 

R«  2.700000EO 
S-  14  DP-REAL 
T-  1.830000E-4 
U-  298 
V-  117 
W-  213 
X-  167 

y-  14  DP-REAL 
2-  14  REAL 


1 

Pr««a  Esc  when  dons  1 

F3  1>  \ 

1, 38E-16*U*F1/ (3*pi*T*2) 

2> 

3> 

4> 

5> 

USB 

PgUp  and  PgDn  to  edit 

other  functions 

next 

*  fUilwC^WilCf 

store  pedlt 

status  print 

dup 

drop 

Bwap 

roll  pi 

Bin 

asln 

alnh 

asinh  inv 

con 

acos 

coah 

acosh  neg 

tan 

atan 

tanh 

atanh  abe 

oxp 

In 

lOexp 

log  sgrt 

B  naxn  nenu  up&xon* 

Help 

Graphics 

Functions 

File 

Proc 

File  I/O 

Polya 

Wave 

Proc 

Curve  Fit 

Diff  Eqs 

Wave 

Gen 

Stats 

save/Exit 

'  Paraaaters 

A-  1000 

B- 

2.5287B1E0 

C- 

0 

D- 

O.OOOOOOE-l 

E- 

O.OOOOOOE-l 

F- 

O.OOOOOOE-l 

G- 

O.OOOOOOE-l 

H- 

O.OOOOOOE-l 

1- 

O.OOOOOOE-l 

Variables 


rho 

R- 

2.700000EO 

beta 

S- 

14  DP-REAL 

eta 

T- 

1.830000E-4 

T 

U- 

298 

LI 

V- 

117 

L2 

w- 

213 

L3 

X- 

167 

data 

y- 

14  DP-REAL 

aizus 

z- 

14  REAL 
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. Function  Definitions  ■  — — 

Press  Esc  when  done 

l>  \  pi*F2 
2> 

3>  \  B*ein(pi/B)*(A'*(F3*(B-l)))‘(l/B) 

4>  / 

5> 

Use  PgUp  and  PgDn  to  adit  other  functions 
—  Calculator  Funotions  — i— — — 


next 

store 

pedit 

status 

print 

due 

drop 

swap 

roll 

pi 

sin 

asln 

alnh 

asinh 

inv 

COB 

acos 

cosh 

aeosh 

nag 

tan 

atan 

tanh 

atanh 

aba 

axp 

In 

10  exp 

log 

sqrt 

Main  Menu  Options 


Kelp 

File  Proc 
Wave  Proc 
Wave  Gan 


Graphics 
File  I/O 
Curve  Fit 
Stats 


Functions 
Polys 
Dlff  Eqs 
Sava/Exlt 


Paraneters 

1000 

2.5237-120 
0«  0 

0-  O.OOOOCOE-1 


E- 

0. OOOOCOE' 

F- 

0. OOOOCOE' 

G» 

0. OOOOCOE' 

H- 

0. OOOOCOE' 

I« 

0.  OOOOCOE' 

i-«  Variables  ^ 

R-  2.70000:i0 

rho 

beta 

s- 

14  DP-REAL 

eta 

X- 

1. 8300002-' 

T 

u- 

298 

LI 

v» 

117 

L2 

w- 

213 

L3 

X- 

167 

data 

Y- 

14  DP-REAL 

alzas 

z« 

14  REAL 

— —  Function  Definitions  .  . — 

Prase  Esc  when  dons 
1>  \  2*F2/F4*0.010fia+F2/ (V*tanh(F4/2) ) 


PgUp  and  PgDn  to  adit  other  functions 
■  Calculator  Functions  i  — ^ 


store 

drop 

asln 

aeos 

atan 

In 


pedit 

swop 

■inh 

cosh 

tanh 

lOsxp 


Main  Kenu  Options 


Help 

Fils  Proc 
Wave  Proc 
Wave  Gan 


Graphics 
Fils  I/O 
Curve  Fit 
Stats 


status  print 

roll  pi 

asinh  inv 

aeosh  nog 

atanh  abs 

log  sqrt 


Functions 
Polys 
Diff  Eqs 
Sava/Exlt 


•  Parameters 

A-  1000 
B-  2.S237S123 
C-  0 

D--  O.OOCOC:2-1 
E-  O.OOOOCCE-1 
F-  O.OOOOC5S-1 
G-  0.0000c:2-l 
H-  O.OOOOCSS-1 
I-  'O.OOOOCtE-l 

>  variables 


rho 

R- 

2.70000020 

bata 

S- 

14  DP-RE.\L 

eta 

T- 

1.8300002-4 

T 

U- 

298 

LI 

v- 

117 

L2 

w- 

213 

L3 

X- 

167 

data 

Y- 

14  DP-REAL 

sizes 

Z- 

14  REAL 

Curve  Fitting  to  an  Arbitrary  Function 


A  -  <  1000  > 

B  -  <  2.52a7BlE0  > 

C  ■  c  0  > 

0  -  <  O.OOOOOOE-1  > 

E  -  <  O.OOOOOOE-1  > 

P  «  <  O.OOOOOOE-1  > 
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Fit  B?  <  y  > 
Fit  C?  <  H  > 
Fit  D7  <  N  > 
Fit  E?  <  N  > 
Fit  F?  <  N  > 
Fit  G?  <  N  > 
Fit  H?  <  N  > 
Fit  I?  <  N  > 
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Fit:  Usar  Function 
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MIE  Sensor  Response  Weighting 


TheMIE,  Inc.  monitoring  instmments,  RAM-1  and  RAS-2  are  very  nonlinear  in 
their  particle  size-dependent  response.  In  order  to  correct  the  measured  mass  response  for 
specific  size  spectra  the  manufacturer-supplied  calbration  curve  must  be  used.  If  a  mono-disperse 
challenge  is  being  used,  then  a  relative  response  may  be  read  directly  from  the  curve.  How'ever, 
we  are  involved  for  the  time  being  with  various  poly-dispersed  challenges  and  size-dependent 
penetration  effects.  Hence  the  nominal  broad-spectrum  response  of  the  instruments  must  be 
corrected  for  the  nonlinear  response. 

A  Mathcad  utility  program,  MIE.MCD,  has  been  written  whose  purpose  is  to 
weigh  a  given  challenge  spectrum  by  the  manufacturer's  response  curve  and  correct  the  nominal 
analog  output  accordingly.  The  program,  shown  in  the  next  two  pages,  allows  two  differen? 
methods  of  presenting  the  challenge  spectrum;  either  by  direct  input  of  the  spectrum  by  size  bins 
or  by  lognormal  distribution  parameter  estimates. 

The  program  first  plots  the  MIE  calibration  distribution  given  as  the  emi  vector 
against  the  particle  sizes,  d;.  Below  the  calibration  plot  is  the  particle  size  spectrum,  given  as 
vector  nj.  Aficr  the  size- weighted  efficiency,  eff,  is  calculated,  the  nominal  millivolt  responses  of 
the  instruments  is  calculated.  These  are  the  values  supplied  by  the  manufacturer.  Their  validity 
must  be  checked  periodically  for  each  instrument  by  gravimetric  calibration  with  well-defined 
particles. 


The  RAS-2  sensor  has  a  nominal  analog  output  of  16  mv/mg/m^  when  calibrated 
to  the  manufacturer's  specifications.  The  RAM-1  monitor  has  three  panel  scale  settings  at  0-2,  0- 
20  and  0-200  mg/m^  with  0-10  volt  output.  In  order  to  prevent  confusion,  the  program  supplies 
all  three  numbers. 

The  option  of  using  the  lognormal  distribution  is  shown  on  the  second  page.  The 
distribution  parameter  estimates  are  initially  user-supplied.  This  program  segment  may  be  used  as 
a  rough  tool  for  manually  optimizing  the  parameters  against  empirical  data  by  comparing  peak 
position,  width,  intercept  and  slope  values.  The  distribution  is  shown  in  the  plot.  Additional 
statistics  based  on  the  lognormal  parameter  estimates  are  available  from  the  Hatch-Choate 
conversion  table  (Reist,  1984).  M  in  the  first  segment,  the  instrument  sensitivities  are  corrected 
by  the  weighted  efficiency  and  displayed  in  the  final  table. 
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MIE  SENSOR  EFFICIENCY  CALIBRATION 
MIE.MCD 


Page  1 


microns 

RAS-2 

Sensor; 

(16mv/mg/M3) 

RAS2  s  eff-16 

RAS2  =17.6 

mv/ing/m3 

RAS-1 

Monitor: 

scale  (nig/m3) 
0-2 

RAS12  »  eff-5000 

RAS12  =  5512.3 

• 

0-20 

RAS120  s  eff'500 

RAS120  =  551.2 

■ 

0-200 

RAS1200  s  eff*  50 

RAS1200  =  55.1 

■ 
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Page  2 


Approximation  for  lognormal  distribution: 


Hatch-Choate  conversions: 


definition  (dp)  p 

Mode  -1 

Geometric  mean  or  mode  0 

Arithmetic  mean  0.5 

Diameter  of  average  area  1 

Diameter  of  average  mass  1.5 

Surface  median  diameter  2 

Surface  mean  diameter  2.5 

Volume  (mass)  median  diameter  3 

Volume  (mass)  mean  diameter  3.5 


Data  identification: 
Talc  nozzle/inside 
for  1221T2 


Distribution  parameters; 
8igOB2.15  meanGB.995 

chosen  p:  p  b  3 

dp  ■  meano-expCp-  In  ( sigG)  0 

dp  =  5.771  microns 


lognormal  distribution  calculation: 


Inorm^ 


V(2'«)  •  In ( SigG) 


•exp 


^„(  )] 

n 

XmeanQ/  1 

v 

\  In ( sigO)  / 

/ 

-2  J 

dj  Inona, 


0. 

053 

0. 

154 

0. 

226 

0 

29 

0. 

409 

0. 

509 

0. 

641 

0. 

696 

0. 

698 

0 

67 

0. 

626 

0. 

574 

0. 

521 

0. 

399 

0. 

301 

0. 

227 

0. 

172 

0. 

101 

0. 

061 

0, 

039 

0. 

025 

0. 

017 

0. 

Oil 

0. 

008 

0. 

006 

0. 

003 

0.1 

0.15 

0.175 

0.2 

0.25 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

1 

1.25 

1.5 

1.75 

2 

2.5 

3 

3.5 

4 

4.5 

5 

5.5 

6 

7 

microns 


eff  ■ 


2  InorTOi'sm^’dj 

i  _ 

lnorm^«d^ 

i 


eff  s  1.388 


RAS2  «  eff-16 

RAS2  = 

22.2 

iiiv/mg/m3 

RAS12  >  eff-5000 

RAS12  = 

6940 

■ 

RAS120  ■  eff-500 

RAS120 

694 

■ 

RAS1200  •  0ff-5O 

RAS1200 

=  69.4 

■ 
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C.6.  Aperture  Flow  Characteristics 

The  air  exchange  rate  of  the  enclosure  during  a  test  is  determined  by  the  flowrate  of 
air  through  the  aperture.  In  some  cases  the  flowrate  may  be  measured  directly  by  any  suitable 
instrument.  If  no  instrument  is  available  for  the  range  involved,  it  may  be  estimated  using  the 
crack  flow  equation; 

Q=CA^-i.dpf 

Use  of  the  equation  requires  estimation  of  the  effective  leakage  area  and  the  exponent 
n.  A  Mathcad  program  has  been  written  to  estimate  these  parameters  from  given  flowrate  and 
pressure  differential  data.  The  utility  program  is  applicable  to  any  test  for  which  the  flow  through 
leaks  must  be  estimated,  whether  a  single  leak  or  a  combination  of  leaks  in  an  unknown 
configuration  such  as  may  be  encountered  with  a  vehicle  or  shelter. 

In  the  program  QPDOC.MCD  on  the  following  page,  the  above  equation  is  applied 
to  finding  the  parameters  for  the  0. 1mm  x  40nim  x  840mm  aperture  used  in  the  verification  tests. 
Since  the  equation  is  more  familiar  with  American  units,  the  constant  value  400S  has  been 
retained  and  the  units  adjusted  for  titers/minute  and  mmWG.  The  natural  logarithms  of  the 
flowrate  and  the  pressure  differential  are  fitted  by  the  linear  least  squares  optimization  function  to 
provide  estimates  of  the  intercept,  m,  and  slope,  n.  The  correlation  coefficient  is  also  calculated 
(corr).  The  effective  area  is  then  calculated  and  given  in  cm^.  Dividing  by  the  length  of  the 
aperture  gives  an  effective  height,  assuming  zero  residual  leakage  area  in  the  enclosure. 

Measurements  of  flowrate  and  differential  pressure  in  the  very  low  range  of  interest 
were  carried  out  using  the  Gillibrator  Primary  Flow  calibrator.  The  blower  was  not  capable  of 
overcoming  resistances  in  the  flow  line  with  any  other  available  flowmeter  so  maesurements  were 
done  very  accurately  in  the  range  0.37-0.57  mmWG  and  results  were  extrapolated  to  -1  mmWG. 

A  slope  of  n=0.5 12  with  a  correlation  coefficient  of  0.992  was  calculated.  The 
effective  area  of  1.24  cm^  indicated  an  effective  aperture  height  of  0. 148  mm.  Under  the 
assumption  of  zero  residual  leakage,  this  difference  from  the  nominal  0. 1  mm  height  could  include 
some  surface  roughness  and  inaccuracy  in  the  simple  feeler  gauge  measurement  which  served  to 
fix  the  aperture  height.  These  optimized  parameters  were  then  used  in  the  protection  factor 
model  to  calculate  the  air  exchange  rate. 
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CALCULATION  OF  LEAKAGE  AREA  CHARACTERISTICS 

QPDOC.MCD 

Characteristic  curve  for  nominal  0.1  mm  x  40  mm  x  840  mm  aperturel215Tl  data 
Number  of  data  points  N  =  5  i  =  0..N  -  1 


data  values: 


Qj  -  Pmmi 


21 

0.55 

19,5 

0.47 

20,5 

0.52 

17,4 

0.37 

22 

0.57 

1pm  mmWG 


calculate  effective  area: 


Convert  pressure  to  in  H20: 
Pnimj 

p  := - i 

‘  25.4 

W,  :=  ln(P,) 

calculate  slope  and  intercept; 
m  :=  intercept(lnP,lnQ)  m“  1.676 

n  :=  slope(lnP,lnQ)  n  “0.512 

corr(lnP,lnQ) -0.992 


convert  flow  to  cfrn: 

q.»A 

‘  28.3 


InQi  :=ln(Qj) 


area  =  exp(ni  -  ln(4005))' 144-6.452  area  -  1.24  sq.  cm. 

calculate  effective  height: 

height  height -0.0148  cm. 

84 

extrapolate  flow  to  -1  mmWG: 

Qex  :=  exp(m)  0.03937”28.3  The  flow  is  Qex  “28.839  1pm  at  ImmWG 


100 


exp^nlnPj  +  mj'28.3 
Q.28.3 


Pmni| 

Pressure  [mmWG] 
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APPENDIX  D 

Piezorcsi stive  Pressure  Transducer  Calibration  Setup 

The  piezoresistive  pressure  transducers  acquired  for  the  experimental  setup  provide 
stable,  linear  voltage  outputs  over  a  wide  range  of  differential  pressures  with  exceptionally  high 
frequency  response.  Their  characteristics  are  presented  on  the  following  two  pages.  With  the 
maximum  amplification  provided  by  the  Endevco  Model  106  Piezoresistive  Conditioner,  the 
tranducers  could  be  used  at  a  sensitivity  of  about  70  mv/mmWG  differential  pressure. 

Since  the  transducers  are  intended  for  use  at  very  low  pressure  differentials  (fi'om  1 
mmWG)  it  is  necessary  to  calibrate  them  against  a  reliable  standard,  The  water-filled  U-tube  can 
be  considered  a  primary  standard,  by  definition.  An  inclined  water  manometer  can  be  used  for  the 
lowest  pressure  range.  In  order  to  readily  calibrate  the  transducers  over  the  nominal  range  of  the 
transducers  (1  psig)  without  recourse  to  a  standards  laborato*~;,  a  simple  setup  was  constructed  as 
shown  in  the  following  schematic  diagram.  The  U-tube  is  about  900  cm.  in  length,  titled  halfway 
with  water.  A  i  ubbor  tube  leads  from  one  arm  of  the  tube  to  a  T-junction.  One  end  of  the 
junction  is  firmly  attached  to  the  threaded  sleeve  supplied  by  the  manufacturer  wl  'ch  seats  onto 
the  body  of  the  transducer  with  a  rubber  0-ring  The  other  arm  of  the  junction  is  connected  to  a 
rubber  bulb.  The  transducer  is  connected  to  the  conditioner  according  to  the  manufacturer's 
instructions  and  allowed  to  warm  up. 

Constant  pressure  on  the  rubber  bulb,  maintained  by  a  laboratory  clamp,  produces  a 
change  in  the  water  level,  measured  with  a  metric  ruler,  and  the  corresponding  change  in  the 
transducer  output,  indicated  on  a  millivoltmeter.  A  calibration  chart  foi  the  range  of  interest  can 
readily  be  constructed  directly  in  terms  of  the  standard  pressure  unit. 


ENDEVCO* 

PRESSURE  TRANSDUCER  TEST  REPORT 


MODEL  8510B-1 

SERIAL#  AAAA7 

Range 

1 

psig 

Sensitivity 

242.9 

mV/psi 

Excitation 

10.00 

Vdc 

Zero  Pressure  Output 

-1 

mV 

Full  Scale  Output 

243 

mV 

Non-Llnearlty 

.51 

%FSO 

Hysteresis 

.07 

%FSO 

Non-Repeatabitity 

.08 

%FSO 

Combined  Lin.,  Hyst.,  &  Rep.* 

.52 

%FSO 

Thermal  Zero  Shift 

.  16 

%FSO 

Zero  Shift  After  3XFSO 

.01% 

3XFSO 

Thermal  Sensitivity  Shift 

.93 

% 

Input  Resistance 

1930.2 

n 

Output  Resistance 

1682. 6 

a 

Isolation  Resistance 

>100 

MO 

•RSS 


NON-LINEARITY 
(%  FSO) 


0  20  40  60  80  100  %  FSO 


+5 


THERMAL 
ZERO  0 
SHIFT 
(%  FSO) 


0  so  100  IBO  200  'F 
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EIMOE  VCO' 

PRESSURE  TRANSDUCER  TEST  REPORT 


MODEL  85108-1 

SERIAL  it  AAAA8 

Range 

1 

psig 

Sensitivity 

3B2.7 

mV/psi 

Excitation 

10.00 

Vdc 

Zero  Pressure  Output 

1 

mV 

Full  Scale  Output 

2B3 

mV 

Non-Linearity 

1.36 

%FSO 

Hysteresis 

.04 

%FSO 

Non-Repeatability 

.OB 

VoFSO 

Combined  Lin..  Hyst.,  &  Rep.* 

1.36 

oFSO 

Thermal  Zero  Shift 

.94 

jFSO 

Zero  Shift  After  axFSO 

0% 

aXFSO 

Thermal  Sensitivity  Shift 

2.76 

% 

Input  Resistance 

1412.6 

0 

Output  Resistance 

1304.6 

n 

Isolation  Resistance 

>100 

MO 

•RS3 


Form  ?6707 


100. 


PRESSURE  TRANSDUCER  CALIBRATION  SETUP 
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APPENDIX  E 

TSI 3302  Diluter  Operation  with  the  LAS-X  Laser  Spectrometer 

In  certain  instances  the  challenge  cloud  concentration  is  liable  to  be  too  great  for  the 
LAS  -X  Laser  Spectrometer  to  handle  without  distorting  the  data.  This  being  so,  a  TSI  3302 
Diluter  with  x20  dilution  has  been  adapted  for  this  purpose.  The  diluter  is  an  accessory  to  the  TSI 
APS  33  Aerodynamic  Particle  Sizer  using  a  S  1pm  total  flowrate  supplied  by  the  sizer.  The  diluter 
works  in  a  closed  system  where  a  small  sample  of  the  aerosol  is  diluted  with  filtered  air  from  the 
input  air  flow.  The  aerosol  flow  is  monitored  by  the  pressure  drop  across  a  laminar  flow  capillary 
tube.  The  dilution  ratio  is  determined  by  a  valve  in  the  clean  air  flow  stream.  The  diluter  is 
normally  placed  on  top  of  the  APS  33  with  the  particle  flow  stream  facing  downwards. 

In  order  to  use  the  diluter  with  the  LAS-X,  an  adaptor  had  to  be  devised  which 
would  feed  the  LAS-X  input  stream  isokineticaily  from  the  diluter  outflow  stream.  The  LAS-X 
samples  at  5  cc/sec  whereas  the  diluter  operates  at  5  1pm.  Hence  the  diluter  output  flow  diameter 
had  to  be  adjusted  to  match  the  flow  velocity  into  the  LAS-X  as  closely  as  possible.  The  diagram 
on  the  next  page  schematically  depicts  the  arrangement  devised  to  overcome  this  problem. 

The  solution  is  based  on  a  1/8”  I  D.  plastic  tube  suspended  within  a  S/16”  I.D.  plastic 
tube.  A  teflon  T-junction  supports  the  inner  tube  and  permits  the  diluted,  reduced  flow  to  be 
sampled  by  the  LAS-X  internal  pump.  The  excess  flow  is  pumped  out  of  the  third  arm  by  the 
main  diaghram  pump  which  is  adjusted  to  maintain  the  proper  How  through  the  diluter.  This 
function  normally  is  fulfilled  by  the  APS  33. 

This  is  obviously  a  temporary  solution,  but  it  is  effective.  Particle  count  reduction 
of  factor  20-22  were  obtained  by  carefully  adjusting  the  diluter  pressure  ratio  according  to  the 
manufacturer's  specifications.  It  should  be  noted  that  the  pressure  ratio  is  very  sensitive  to  the  air 
sampling  geometry  at  the  inlet  of  the  diluter.  If  a  long  sampling  tube  is  used,  the  pressure  ratio 
must  be  adjusted  carefully  to  the  new  situation. 

As  noted  in  the  text,  the  penetration  efficiency  as  a  function  of  particle  size  of  the 
diluter  capillary  is  given  as  a  graph  in  the  manufacturer's  instruction  manual.  The  graph  is 
normally  stored  in  the  APS  software  for  automatic  use.  In  order  to  make  use  of  the  correction 
curve  in  our  utility  programs,  the  curve  was  fitted  to  a  polynomial  function  with  ASYSTANT. 

The  resulting  best  fit  parameters  are  included  in  the  LAS-X  data  analysis  program.  The  best  fit 
equation  for  the  x20  diluter  is  as  follows: 

£  =  1.06-0. 024-D-0.002-D* 

where  E  is  the  penetration  efficiency  and  D  is  the  particle  diameter  in  microns.  Penetration 
degradation  begins  at  approximately  2  microns  and  reaches  60%  efficiency  at  10  microns. 
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ADAPTATION  OF  TSI  3302  DILUTER  TO  LAS-X  SAMPLING 
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Dissemination  System: 
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Figure  F.4.  Sonic  Nozzle  Disseminating  Talc 


Figure 


Figure  F.6.  Test:  Enclosure  View  from  Wifhin  Facility 
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F.8.  Test  Instminentation  Inside  Enclosure 


